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Abstract 



i-^ I The rigorous study of spectral stability for strong detonations was begun by J. J. Erpenbeck in 

^p^i [Erl]. Working with the Zeldovitch-von Neumann-Doring (ZND) model, which assumes a finite 

C"! ' reaction rate but ignores effects like viscosity corresponding to second order derivatives, he used 

. a normal mode analysis to define a stability function V{t, e) whose zeros in 3?t > correspond to 

multidimensional perturbations of a steady detonation profile that grow exponentially in time. 
Later in a remarkable paper [Er3] he provided strong evidence, by a combination of formal and 
rigorous arguments, that for certain classes of steady ZND profiles, unstable zeros of V exist for 
perturbations of sufficiently large transverse wavenumber e, even when the von Neumann shock, 
regarded as a gas dynamical shock, is uniformly stable in the sense defined (nearly twenty years 
. later) by Majda. In spite of a great deal of later numerical work devoted to computing the 

If^ ' zeros of V(t, e), the paper [Er3j remains the only work we know of that presents a detailed and 

. convincing theoretical argument for detecting them. 

' The analysis in [Er3] points the way toward, but does not constitute, a mathematical proof 

. that such unstable zeros exist. In this paper we identify the mathematical issues left unresolved 

in [Er3] and provide proofs, together with certain simplifications and extensions, of the main 
conclusions about stability and instability of detonations contained in that paper. The main 
mathematical problem, and our principal focus here, is to determine the precise asymptotic 
behavior as e — > cx) of solutions to a linear system of ODEs in x, depending on e and a complex 
frequency r as parameters, with turning points on the half-line [0,oo). 
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1 Introduction 

The most commonly studied model of combustion is the Zeldovitch-von Neumann-Doring (ZND) 
system ()A.14p , which couples the compressible Euler equations for a reacting gas (in which pressure 
and internal energy are allowed to depend on the mass fraction A of reactant) to a reaction equation 
that governs the finite rate at which A changes^ In three space dimensions with coordinates (x, y, z) 
a steady planar strong detonation profile is a weak solution of this system depending only on x 
with a jump (the stationary von Neumann shock) at x = 0. Without loss of generality we study 
profiles of the form w = {v,u,0,0, S, X), where v is specific volume, n > is the x-component of 
particle velocity, and S is entropy. The solution is constant and supersonic {u > cq, where cq is the 
sound speed at x) in x < 0, the quiescent zone, and satisfies a nonlinear system of ODEs in the 
subsonic reaction zone x > 0. In order to be a weak solution in a neighborhood of x = it must 
satisfy an appropriate Rankine-Hugoniot condition at x = 0. There is a well-defined limiting state 
Woo = lima;_>.oo w{x), which represents a state of chemical equilibrium, and the range of u on [0, oo) 
is a compact subinterval of x > 0. 

If one perturbs this solution, say at time t = 0, with a multidimensional perturbation w{x, y, z), 
the profile w{x) and the front defining the von Neumann shock will change with time, their de- 
scription at time t now being given by functions W{t, x, y, z) and x = ip{t, y, z) satisfying the ZND 
system, a Rankine-Hugoniot condition, and the initial conditions 

(1.1) W\t=o = w{x) + w{x, y, z), i;\t=o = 0. 

To a linear approximation the evolution of the perturbation is governed by the linearization of the 
ZND system with respect to both W and about the steady profile w{x) and steady front x = oH 
General perturbations can be represented by Laplace-Fourier transform in (t, y, z) as superpositions 
of oscillatory ones, so it is natural to look for normal mode solutions of the linearized problem of 

^Without changing the analysis below one can allow n reactants, in which case A is an 7i-vector; for convenience 
we take n — 1. 

^In [Er3| and in his work on non-reactive shocks [Er4| . Erpenbeck appears to have been the first to recognize the 
importance of linearizing with respect to both the state W and the front ^. 
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the form 



(1.2) W{t, X, y, z) = e*^+^°^+*^^(/.(x, r, a, (3), ipit, y, z) = e^^+iay+i^z^^^^ 

where (a, /3) G and r G C with 5Rt > 0. Solutions (|1.2p with Sir > and decaying as x — )• ±00 
correspond to perturbations that grow exponentially in time. The spectral stability problem is to 
determine whether such frequencies (r, a, (3) exist, and if so, to locate them as well as possible. 

The rigorous study of the spectral stability problem for ZND detonations was begun by Er- 
penbeck in [Erlj . Instead of working with the 6x6 system of equations satisfied by the functions 
(j) = {v,u,Uy,Uz, S, X) and ^ in (jl.2p . he studied the reduced 5x5 system satisfied by ^ and the 
new unknown 

(1.3) ^■.= iv,u,U,S,X) 
whera3 



(1.4) U{x,T,a, (3) := {aiiy + Puz)/e, and e = a/o^ + /J^. 

and the other components of (p are exactly as before. 

The coefficients of the reduced system depend only on {x, r, e) so, dropping the tilde on (p, 
Erpenbeck writes the unknowns as (p{x,T,e) and ^(r, e). The 5x5 system of ODEs on 2; > 
satisfied by the functions (p and ^, together with the linearized jump condition, are given in formula 
(4.1) of [Erlj . The equation on x > is 

(1.5) ^ = P(x,T,e)0 + /(x,r,6) 

ax 

where the exact form of / need not be specified now and 

(1.6) P{x, T, e) = -A-\x)[tI + ieAyix) + B{x)], 

for matrices A^, Ay and B given in (|A.2ip . The x-dependence of these matrices enters entirely 
through the profile w{x). 

In [Erl] Erpenbeck defined a stability function V{t, e) ()A.15P whose zeros in the right half 
plane 5Rr > ("unstable zeros") correspond to perturbations of the steady profile w;(x) that grow 
exponentially in time. Later in a remarkable paper [Er3] he provided strong evidence, by a combi- 
nation of formal and rigorous arguments, that for certain classes of steady ZND profiles, unstable 
zeros of V exist for perturbations of sufficiently large transverse wavenumber e, even when the von 
Neumann shock, regarded as a gas dynamical shock, is uniformly stable. Much numerical work has 
been devoted to locating zeros of V; see, for example, |Er21 IBZ21 \LS\ IShl ISSj IHuZj and references 
therein. The argument of |Er3] . though mathematically incomplete and at times incorrect in ways 
that we describe below, is the only convincing theoretical argument we know of that actually de- 
tects and locates unstable zeros. Our goal here is to provide a mathematically rigorous proof of 
the instability (and stability) results of [Er3J . The main mathematical problem is to determine the 
precise asymptotic behavior as e —t- 00 of solutions to a linear system of ODEs in x depending on 
e and the complex frequency r as parameters, with turning points x* on the half-line [0,oo). 

The computation of V requires the evaluation within the reaction zone of the solution 9{x, r, e) 
of the homogeneous transposed equation 

(1.7) — = -P\x, T, e)9 on X > 



Thus, e/27r is the transverse wavenumber. 
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which satisfies the condition that 6 remains bounded for fixed (r, e) with > as x — )• oo. As we 

will see, this condition determines 9 uniquely up to a constant multiple; morever, for 5Rr > the 
solution 9 decays exponentially to zero as x — )■ oo. 
As in (Er3] we decompose r as 

(1.8) T = Ce + i^, 

where C £ {-^ S C : > 0}, S {z S C : > 0, \z\ < R}, and e > is large. In much of what 
follows z/ will play no essential role and one can take R = 0. However, in cases where unstable zeros 
of V{t, e) do exist, the freedom to vary the parameter is useful in proving the existence of those 
zeros, and for that one must take R sufficiently largeQ Thus we can rewrite equation (|1.7p as 

d9 

— = (e^Q + ^A9 where 
ax 

(1-9) <^o{x,0 = {A-Hx)■{CI + ^Ay{x))Y 

z.) = {A-\x) ■ viy + {A-\x)B{x)Y. 

The eigenvalues of the 5x5 matrix $0 (|A.23p play a crucial role in all that follows. They are 

(1.10) /ii = -K« + s)/7/n, H2 = -k{kC - s)/r]u, = = = (/u, 

where with Cq = —v'^py{v, S, A) 



(1.11) s(x,C) = ^C + clr], k{x) = Y^l - 7? = u/cQ. 

Here the square root defining s, regarded as a function of is taken to be the positive branch with 
branch cut the segment [—icQ^Jri,icQ^J^if\ on the imaginary axis. Thus, in particular, we have 

s = \s\ when + CqT/ > 

(1.12) s = i\s\ when + CqI? < and ( = i\C,\ 

s = —i\s\ when + c^rj < and C = — ^ICI- 

One checks that only ^1 has, for ?R.( > 0, negative real part. The corresponding eigenvectors are 
the respective columns of the matrix 

/ 

(1.13) T(x,C) = 

\ 

where m = - is the mass flux. 

V 

On any subinterval I of [0, 00) where (for fixed C) the transformation T(x, () is invertible, we 
set 9 = T{x, Qtt and obtain the system with diagonal leading term 

^ = {eD + E)tt, where 

(1.14) 

D :=diag(/ii,...,/U5) and E{xX,i^) ■.= T~^^iT -T'^ — . 

ax 
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^Since e = |(Qf,/3)|, one must allow C, to vary over all of 5R(" > to cover all possible triples (r, a,/3) under 
consideration. The minimum size of R necessary for detecting zeros of V is given in part (d) of Theorem 15.31 
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On such subintervals one construct^ approximate solutions to order e of (|1.14p associated to 
each of the eigenvalues /Uj of the form 

(1.15) 7r,(x, T, e) = e^'^»(-'f)+^'(-'f''^)[/,o(x, C, i^) + e'^faix, C,J^) + --- + C, ^)], 
where 

(1.16) h^{x,C) = [ l^.{x',C)dx' 



JO 

and fifi is the unit vector with i-th component 1. This is done by substituting the expansion 
into (jl.l7p and solving the equations obtained by equating coefficients of equal powers of e. For 
z = 1, 2 we will use the formulas 

(1.17) kiix,C,'^)= / Eii{x\Q,u)dx' 

where Cyu is a contour that lies on the real axis, except for short excursions into the upper or 
lower half plane to avoid singularities of En ()A.24p . The choice of depends on the profile w 
and is explained in sections 14.11 and 14.21 We do not need explicit formulas for the other quantities 
appearing in (jl.lSp for i = 1, . . . , 5. Sufficient conditions for approximate solutions of this type 
to be 0{e~^"^^^^) close in relative error to true exact solutions of (|1.14p for e large are given in 
Theorem 13.11 

Prom the formulas (jl.lOp we see that for any fixed value of the eigenvalues and fj,2 are 
distinct except at x values where s"^ = ^'^ + CqT] = 0; at such values the first and second columns 
of T are parallel. The eigenvalues ^2 and fis are distinct except at x values where C, = and then 
the second and third rows of T are clearly parallel. For all other values of x the matrix T(x, Q is 
invertible. The special values x* = x^:{C,) where T is singular are referred to as "turning points". 
These points determine a family of subintervals X C [0, 00), and on each of these subintervals the 
approximate solutions described above can be constructedjj 

A complex number C, with IkQ > is defined in |Er3j to be of Class III or Class II respectively, 
when there exists G [0,oo] such that s(x^,,(^) = or C = n(x*)0 All other C, are said to be of 
Class I. Thus we have 

Class III = {C : = and min(cor/2) < |C| < max(co772)} 

X X 

(1.18) Class II = {C : 9^ = and minu < C ^ maxn} 

X X 

Class I = {all remaining ^ S C with ?R.(^ > 0}. 

Class III (resp. II) consists of two (resp. one) bounded closed interval(s), and the minima appearing 
in (jl.lSp are positive. 

Notations 1.1. When working with Class III values of Q we will usually suppose C, = i\C,\- Similar 
results hold with similar proofs when C, = —i\C,\, hut certain formulas change slightly. Thus, for 
some statements it is helpful to define 

(1.19) Class III+ = Class Illn {C = i\C\}. 



^See for example Chapters 5 and 6 of Coddington and Levinson ICLI. 

®The family of subintervals depends on but we normally suppress ^ in denoting both the intervals and the a;*. 
''When Xt = 00, the obvious limits are intended here. For a given £ Class III or Class II, there may be more 
than one x, such that 5(2;*, Q = ot ( = M(a;«). 
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A simpler form (j5.3p of the stability function V{t, e), based on the choice of a "good unknown" 
that removes the front from the linearized ZND system in the reaction zone, was derived in |C JLWj . 
That form involves no integrals bj as in ()A.16p and was shown in [CJLWj to be equal to Erpenbeck's 
stability functionlfl The formula (|5.3p for V{t, e) depends on the value at a; = of the exact decaying 
solution 9{x,T,e) of (|1.7p . Our assumptions guarantee that 6 and P*(x,t, e) are real analytic in 
X on [0,00), so for fixed (r, e) the function 6 can be extended analytically as a solution of (|1.7p 
for z in some complex neighborhood of [0, 00). For a given C, let us suppose that there are a finite 
number of turning points {N = for C, G Class I): 

(1.20) < xn* < xn-1* < • • • < xu, 

The strategy of [Er3] for determining ^(0, r, e) is first to determine 9 on the interval [xi^, + 6, 00) for 
some 5 > 0, then to analytically continue that solution in C \ xu around the turning point xi* and 
up the real axis to X2* + then to analytically continue around X2* in C \ X2*, and to repeat this 
process until x = is reached. It is impossible to determine explicit formulas for the exact solution 
9, so instead one tries to determine explicit formulas that closely approximate the exact solution for 
e large on a certain collection of subsets of complex plane, subsets whose union covers [0, 00) \ T, 
where 7" is a union of small intervals each centered at a turning point. The approximating formula 
will generally change from subset to subset, so there arises a matching problem in passing from a 
given subset to an adjacent one. 

On bounded subsets of the complex plane we use approximating integral formulas of type (jl.lSp 
and linear combinations thereof. Sufficient conditions for such formulas to closely approximate 
exact solutions of (jl.7p (for e large) together with precise error estimates are given in Theorem 13.11 
We will use the notation 9j to denote an exact solution defined on some bounded region Oj C C 
which lies close to an approximate solution 9j = Tttj defined by a formula of type (jl.lSp on a 
possibly larger region Oj. Theorem 13.11 does not apply to unbounded regions; in order to obtain 
an accurate approximating formula on the unbounded interval [xi* + 6, 00) we use a completely 
different approach based on the variable-coefficient gap lemma. Lemma I A. H introduced recently in 
[Z2] . 

Following [Er3j . when considering Class III values of C we focus on profiles w{x) for which the 
quantity CqIj^x) = Cq(x) — n^(x) either increases monotonically for x S [0, 00) (case I), decreases 
monotonically (case D), or increases to a maximum at xm and then decreases (case M). In case M 
our analysis omits treatment of the special Class III+ value of given by 

(1.21) Cm = icorj^ixM)- 

The turning point argument for such Cm in |Er3j is much more involved than that for the other 
Class III values of C; we do not discuss that case here. The turning point arguments for other 
Class III values of C in case M, given in the proof of Proposition 14.81 are a combination of the 
arguments for cases I and D. Erpenbeck shows on p. 1303 of |Er3) that for "A — )• B reactions" (one 
reaction detonations with no back reaction) with an Arrhenius rate law (jl.30p . the cases I, D, and 
M actually occur, the distinction depending on values of the heat capacity ratio, heat of reaction, 
and Mach number of the detonation. His results indicate that for A ^ B reactions the case of a 
single interior minimum does not actually occur. 

The arguments of |Er3j are incomplete in two ways, namely, in the treatment of the unbounded 
interval [xi^, + 5, 00) and in the problem of matching solutions on adjacent regions. To treat the 

* It was also shown to be a nonzero multiple of the ZND Evans function defined in [JLWj . 
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unbounded interval |Er3| considers the limiting problem 



(1.22) _ = _P*(oo,T,e)0 

obtained by evaluating the coefficients of (jl.7p at the values given by limaj^oo '"^(a^), and refers to 
an argument of [Erlj showing that for large x, solutions of (|1.7p lie close to solutions of (jl.22p . 
However, |Erl] is not concerned with behavior as e — )> oo, and this argument does not address the 
possibility that as e increases solutions of ()1.7p remain close to solutions of the limiting problem 
(|1.22p only on intervals [x(e), oo) such that lim^^oo x{e) = oo. Too rapid growth of x(e), for example 
0{e'-"''), would make it impossible to carry out the matching arguments for the remaining interval. 

To obtain an approximating formula for the decaying solution on the unbounded interval, |Er3) 
observes that the integrals defining hi{x,C) and ki{x,<^,i') "become linear in x near the region 
of chemical equilibrium", i.e., as x — )• oo. Since for most choices of (, the eigenvalue ^i(x,C), 
the corresponding eigenvector of ^o{xX), and £'11(2;, C,i^) approach their counterparts, which are 
constant in x, in the perturbation expansion for e large of the eigenvalues and eigenvectors of 
e^>o(oo, C) + ^1(00, v), \Ei3\ concludes that "the e — t- 00 limit and the x — )■ 00 limit are interchange- 
able" , and that on the unbounded interval the exact decaying solution is well-approximated by the 
approximate solutioii^l 

(1.24) ei{x,T,e) = T{xX)ni{x,T,e). 

He then uses 9i in the matching arguments on the remaining finite interval. 

This argument is suggestive but nonrigorous, and it actually leads to an incorrect conclusion. 
In Theorem 12. II we use an argument based on multiple conjugations and the variable coefficient gap 
lemma. Lemma lA.ll to construct the exact decaying solution 6 on the unbounded interval. The 
construction provides at the same time an accurate approximating formula and an error estimate 
(see (2.9)). Comparison of this formula with 61 shows that the solution 9{x,T,e) differs from 
6i{x, T, e) by a complex factor (namely e^"^^''^''^^ in (j2.35p ) that depends on x and grows exponentially 
with X for certain choices of (Ci^^)- Thus, it is not true that 9 ^ 9i (see the preceding footnote); 
in fact 9 ~ M(t, e)0i is not true on any open x-interval for any multiple M(r, e) of 9i. Even for 
choices of (Ci'^) for which the factor e^'^^''^''^^ remains bounded as x — )• 00, the presence of this 
factor means that the matching argument for the unbounded interval and its adjacent bounded 
region has to be done with a carefully chosen multiple of 9i (or alternatively, of 9) and not 9i itself. 

In the matching arguments one uses bases of exact solutions {^1, . . . , ^5} of ()1.7p whose asymp- 
totic behavior as e — )• 00 is known in regions of the complex plane that include portions of the 
adjacent intervals under consideration. In these arguments one must take into account the growth 
rate with respect to e not only of the solutions 9j (to know, for example, which solutions are "dom- 
inant" or "recessive" in certain regions), but also the growth rate of coefficients. For example, if 
9 denotes the exact decaying solution on the interval [xi* + 6, 00) constructed by Lemma lA.H the 
first step in determining the analytic continuation of this solution around xi^ is to expand 9 in such 
a basis: 

(1.25) 9{z, r, e) = ci(r, e)9i{z, r, e) + • • • + C5(r, 6)^5(2, r, e). 

^Erpenbeck writes "6 ~ ^1" for x > xi,. According to the standard definition in asymptotic ODE theory, given 
in [CLj for example, 

(J 

(1.23) 9 ~ 9i ior X G S there exist positive constants C,R such that |& — Si| < — |^i| for a; G S, e > R. 
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The coefficients in (jl.25p are independent of z but depend on (r, e), and a necessary part of a 
rigorous matching argument, a part that is missing in |Er3j . is to keep close track of the growth 
rates of these coefficients with e as one passes from one basis to another. This type of analysis is 
carried out in the proofs of Propositions 14.31 and 14.61 

In case D we provide some arguments that seem to be missing in the analysis of |Er3j . As in 
case I, where there is also just a single turning point x* for class III values of C, after the decaying 
solution Q has been analytically continued just a little to the left of a;^,, there remains the question 
of determining the asymptotic behavior as e — >• oo of the solution on the remaining part of the 
interval [0,x=k — S\. That determination requires knowledge of the asymptotics of all the exact 
solutions in the basis used to express d (I0ni . (071) . The approach of [Er3] is to use the method of 
the parameter problem, Theorem 13. 11 for this purpose. That approach works only in case I, where 
= on [0, — (5]; in case D it turns out to be impossible to use Theorem 13. II to determine the 
asymptotic behavior of all but one of the basis solutions on [0,x* — S\. The reason, explained in 
step 6 of the proof of Proposition 14.61 is that certain differences ?R.{fii — fij) have an unfavorable sign 
on [0, X,, — S]. In step 6 an energy-type estimate is provided along with an analysis of coefficients 
to complete the treatment of this case. Such arguments are needed as well in case M, where there 
may be two turning points associated to a given ( G Class III, for understanding the asymptotic 
behavior of the exact solution on the full interval between the turning points. 

To treat Class III values of ^ that correspond to three or more turning points, matching argu- 
ments beyond those given in this paper may be needed, particularly for the continuation around 
Xfc* when is of increasing type. 

1.1 Assumptions 

Assumption 1.2. The thermodynamic functions appearing in the ZND system ()A.14j) . p (pres- 
sure), T (temperature), AF (free energy increment), andr (reaction rate) are real analytic functions 
of their arguments {v, S, A). 

Assumption 1.3. The steady strong detonation profile w{x) = {v,u,0,0, S, X) is a real-analytic 
function of x in the subsonic reaction zone [0, oo). There exist constants Ci and C2 such that 

II 

(1.26) < Ci < K = — < C2 < 1 and < Ci < li < C2 < 1 /or a// X e [0, 00). 

Co 

Moreover, there exist positive constants C , f3 such that 

(1.27) |'u;(x) - w;(oo)| < Ce"^^. 

We mainly study profiles of type I, D, or M, and will call attention to results that hold without this 
restriction^^ 

Assumption 1.4. The rate function satisfies 

(1.28) rU=o = 0, rx< 0, r,U=o = 0, rs|A=o = 0. 
This assumption is satisfied, for example, by rate functions of the form 

(1.29) r = -kp(t){T)\ 

^''For C, of Class III there is a single turning point a;* in cases I and D. In case M, depending on the choice of 
£ Class III, there can be one or two turning points. 
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where p is density and /c > is a reaction rate constant, such as the Arrhenius rate law 

(1.30) r = —kXexp{—E/ RT) (E' is activation energy) 
considered in the section of [Er3] on A ^ B detonations. 

Remark 1.5. Assumption 11.41 can be weakened to allow for rate functions (such as those incor- 
porating back reactions B ^ A) for which one may not have r„ = and rs = when r = 
(equilibrium) . Inspection of the proof of Theorem 12.11 shows that it suffices there to have for all 
X > 0: 

(1.31) ^{n* - fil){x,C,'^,h) > -6^h + 0{he~'^''), j = 2, . . . , 5 where 5* < /3, 

where /3 is as in (jl.27p . h = 1/e, and the //^ are the eigenvalues of ^o{x,() + /i$i(x,j/). For rate 
functions satisfying ()1.28p we check below that (jl.3ip is satisfied with 5* = 0. It is hard to check 
(jl.3ip directly, but a readily verifiable condition that implies (jl.3ip is 

(1.32) ^{n* - fil){ooX,J^,h) > -6,hj = 2,..., 5 where 5* < /3, 

the corresponding condition for the constant coefficient limiting system $o(c«, C) + ^$i(oo, ly). For 
such rate functions the remaining arguments of this paper are also valid (see also Remark 12. 6p . 

Assumption 1.6. The stability function for the von Neumann shock, defined originally in JEr^]^ 
has no zeros in ?R.( > 0. 

This means that the equation of state of the unreacted explosive is such that the von Neumann 
step-shock would be stable if the reactions behind it were somehow suppressed. This assumption 
is made in |Er3| to allow us to concentrate on instability which arises solely from the reactions. It 
always holds, for example, for step-shocks in ideal polytropic gases. 



1.2 Main results 

The main results of this paper provide a rigorous determination of 0(0, r,e) (and, in fact, of 
6{x, T, e) for x away from turning points) to arbitrarily high accuracy for e large and are contained 
in Theorem 12.11 Proposition 14. 3| and Proposition 14.61 The behavior at x = for e large may be 
described informally as follows. Let ti(C) and t2(C) be the first two columns of the matrix T{0,Q 
prT3|) . With T = Ce + u we have 

For C ^ Class III, 0(0, r, e) = ti(C) + 0(l/e) 

For C G Class III+ \ E, 0(0, T, e) = ti(C) + a(e, C, z^)t2(C) + 0(1A), 

where a is given by (|4.30p and the exceptional set E = {iCo??^/^(0+), ico?/^/^(oo)} in cases I and D. 
In case M, E also contains ico?7^/^(xM)) where xm is the location of the maximum. 

The main instability result of |Er3| . restated here as Theorem 15.31 follows then by an argument 
using Rouche's Theorem and the formula ()5.3p for V{T,e) in terms of 0(0, r,e). In cases I and M 
this theorem provides an explicit sufficient condition (|5.13p . expressed in terms of functions of the 
steady flow variables, for detecting the presence of zeros of V{t, e) in > when e is sufficiently 
large. The unstable zeros occur for a subinterval of Class III values of at certain values of 
with > 0. In case D the result allows one to conclude for all C except the two particular 

^^This stability function turns out to agree with the function -£/i(C) (|A.19|) . described in section [SI and to be a 
nonvanishing muhiple of the Majda stabihty determinant for shocks defined in |M] 
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Class III+ values in (jl.35p that for e > e(C) sufficiently large, V{C€ + v, e) does not have unstable 
zeros. Moreover, for the A ^ B reactions described above, [Er3| shows that for certain values of 
the physical parameters, unstable zeros are actually present in cases I and M. 

Our use of Lemma lA.ll permits certain simplifications and extensions relative to |Er3j . The 
application of the lemma to this problem depends only on having good separation between /xi and 
the remaining set of eigenvalues. Thus, the crossing of fj,2 and /X3 that occurs for Class II values 
of C at the associated turning points are irrelevant for us. In our analysis it is only necessary to 
distinguish two sets of C, namely, Class III and its complement in > 0. In |Er3j the unbounded 
interval could be taken to be [0, 00) only for Class I values of C; here we may take it to be [0, 00) 
for both Class I and Class II. The special Class II values given by 

(1.34) C = u{0+) and C = hm u{x) := u(oo), 

which correspond to turning points "at and at 00", had to be excluded from the analysis of [Er3j , 
but they can be treated in our approach like any other Class II values. Thus, for all non-Class III 
values, 9{0X,u,€) can now be determined in a single step, with no matching required. However, 
like |Er3j we must exclude from consideration the exceptional Class 111+ values 

(1.35) {icor?5(0+),ico7?5(oo)} 

in the matching arguments of section 

Lemma lA.ll also permits us to keep track of how the constants that appear in our estimates 
depend on For example, our proof of Theorem 15.31 shows that for in any compact subset K 
of {{'Sic ^ 0} \ Class III) there exists a positive constant e{K) such that V{T,e) 7^ for e > e{K). 
Obtaining uniform wavenumber cutoffs €{K) has an obvious importance for numerical investigations 
of the presence of unstable zeros. There was no attempt to get uniform cutoffs in [Er3j, and in fact 
such cutoffs cannot be extracted from the arguments given there for reasons already indicated (for 
example, incorrect treatment of the unbounded interval and the need to use matching arguments 
for Class II values of ( due to crossing of //2 and fi^). Moreover, our proof of the cutoff result just 
stated holds even for profiles that are not of type I, D, or M; such profiles are not studied in [El^. 
We also prove statements about uniform choices of e in the parts of Theorem 15.31 that pertain to 
Class III values and instability. 

We introduce a major simplification in the proof of Theorem 15.31 by working with the simpler 
form ()5.3p of V{T,e), instead of with Erpenbeck's original form ()5.ip . The analysis in |Er3 j of the 
asymptotic behavior of the integrals bi and 62 (|A.16p . which is complicated since the integrands 
involve 9, is thus no longer needed. The use of the simpler form of V also gives more rapid 
convergence of the approximate stability function La to L (see (j5.4p : here V = eL), which in turn 
yields wavenumber cutoffs e{K) that are better (i.e., smaller) than those derivable using the original 
form of V. 

1.3 Discussion and open problems 

A rigorous approach to some aspects of the high frequency behavior of the ZND stability function 
y(r, e) was given in |C JLW] . That paper addressed the question of the behavior of the exact 
decaying solution 6 of (jl.7p on the unbounded interval. A different type of gap lemma was used 
there to show that 9{x,T,e) (suitably normalized) is well-approximated by the solution 6L(x,T,e) 
of the limiting problem (jl.22p on unbounded subintervals [x(e),oo), where x(e) converges to 00, 

Although |Er3) does not explicitly exclude icor]i(0+), we believe that his analysis does not cover this case. 
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but not too fast to rule out matching arguments. No comparison was made there between 6l and 
01- 

The present treatment based on Lemma [A . 1 1 provides several improvements over |C JLWj . First, 
we obtain an approximating formula for 6 that is valid on unbounded intervals whose left endpoint 
is independent of e. In Theorem 5.1 of jCJLW) a separate argument using the Tracking Lemma 
of [Z3] was used to prove the absence of unstable zeros for e sufficiently large only for a certain 
proper subclass of non-Class III values (mostly away from 3ft(^ = 0). Turning points, the behavior 
of y(r, e) for Class III values of C, and instability were not treated in [CJLWj. 

Our study completes the mathematical analysis of |Er3| . yielding rigorous general criteria for 
instability of ZND detonations. Despite extensive numerical and asymptotic studies of instability 
in the detonation literature, this represents to our knowledge the only analytic proof of instability. 

Another interesting and closely related problem is to establish conditions for stability of deto- 
nation waves in the small heat-release limit AF — )• 0, in which equations ()A.14p formally decouple, 
and the associated high overdrive limit^^ as carried out in the one-dimensional case in |Z2j . As 
discussed in [Z2], both of these problems reduce to establishing uniform high-frequency stability 
estimates similar to those obtained in the present paper; as noted in jEr2] , the latter problem is 
quite subtle in multi-dimensions. Such a result would, together with the instability results devel- 
oped in |Er3j and this paper, complete the picture of behavior that has emerged from numerical 
lore of transition from stability to instability as heat release or overdrive is varied. 

The study of positive high-frequency stability criteria is further motivated by the needs of nu- 
merical stability analyses, where, to obtain definitive results of stability, it is necessary to perform 
some such asymptotic analysis truncating the computational domain to a bounded set of frequen- 
cies. See |BZH IBZ2t IZ2] for treatments in the one-dimensional case. 

What is needed beyond the analysis of the present paper to obtain positive high-frequency 
stability criteria is to treat the Class III frequencies neglected here and in [ErS] for which turning 
points X* occur at x = or x = +oo, and (for uniform estimates) frequencies in their vicinity, as 
well as the additional frequencies neglected here in case M. The latter should be treatable as in 
|Er3| . We expect that the problem of obtaining uniform estimates as |C| — oc0 can be treated as 
a perturbation of the one-dimensional case using the arguments of [Z2j . A remaining question is 
to obtain estimates that are uniform as C approaches the imaginary axis near Class III values. At 
the moment we have such estimates near Class I values on the imaginary axis. We plan to address 
these questions in a separate work. 

Though they were not needed in order to obtain useful instability results, we regard the treat- 
ment of these remaining cases, and subsequent applications to positive stability results, as a very 
interesting open problem. Another interesting open problem is to explore the implications of these 
and the present analyses for the associated viscous problem, as in \Z1\ ITZ] , 

Remark 1.7. One can try to study the behavior near turning points using oscillatory integrals. 
Near turning points of increasing or decreasing type these have the form associated with fold 
caustics and lead to Airy functions (see also Remark 14. 9p . In future work we will explore whether 
such an approach might yield a simpler, or more efficient, or more unified treatment. 

^■^ As described in [Zll IZ2] . this amounts to a simultaneous large shock-strength and zero heat release/activation 
energy limit; for further description, see [Z2| and Appendix C, [Zl] . 
^^Recall that we prove uniform estimates for bounded |(^| only. 
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2 Exact decaying solutions on [a, oo), a > 

In this section we prove Theorem 12. II for the unbounded interval. The theorem gives an accurate 
approximating formula when e is large for the exact decaying solution 6{x, ^, i/) of (|1.7|) on the 
interval [x^, + S, oo) for Class III values of ( and on the interval [0, oo) for all other values of C in 
3f?C > 0. 



2.1 Construction of exact solutions 

Recall from (fLQl) that 



M^x) = {A^Hx)-{ci + iAy{x))y 

'^iix,u) = {A-\x) ■ i^iy + {A-\x)Bix)Y, 
where the coefficient matrices are defined in section IA.4I So with h = ^ we can write 

(2.2) ^o{x, C) + h<^i{x, v) = ^o{x, C + H + h{Ax\x)B{x)y. 

The eigenvalues of ^o{x, ( + hv) are fj.j{x, ( + hu), j = 1, . . . ,5, for fij as defined in (jl.lOp . 

In preparation for the theorem, we need to examine the effect of h{A~^[x)B[x)Y on the eigen- 
values of ^o{x,C) + h^i{x,i'). Direct computation and the use of Assumptions 11.31 and 11.41 shows 
thaS 

(2.3) {A-\x)B{x)Y = Oie-^"^) + ( ° . ) , where row 5 = (*, *, *, *, -rx/u), 

Vrow o J 

a matrix whose only nonzero row at equilibrium is the fifth. This gives 

(2.4) <^o{xX) + h<^i{x,u) = <^>o{xX + hu) + h( ^ }\ + 0{he-^''). 

Vrow / 

Prom (IA.23P we see that ^q{x,C) is a matrix of the form 

where the eigenvalues of A are /xi, /U2, //a = /X4 = C/""- Thus, the eigenvalues of ^o{x, () + h^i{x, v) 
are 

^l* = ^lJ{xX + hv) + 0{he-^''), j = 1,2,3,4 

(2-6) 

^lt = H3{x, C + hu)-h— + 0(/ie-^^), where rx < 0. 

u 

Here we are using (12.41) and the fact that the eigenvalues of $o + f ^ ^ I : which are semisimple, 

yrow 5 J 

undergo a perturbation of size 0{he~^^) in response to a matrix perturbation of the same size. 



^^The computation is done on pages 114-117 of [Er3| . There is a sign error in the expression for is ■ Wii, which 
should be multiphed by —1, in (A. 15) on pll7 of [Er3] . This correction yields the expression for row 5 in 1)2. 3|) . 



12 



Theorem 2.1. With h = 1/e consider the system (jl.7p 

(2.7) 0' = ^[<^o{oo,C) + h^i{x,iy)]9 
on an interval [a,oo), a>0, and for values of such that 

(2.8) 1^1 (x, C + hu) — fij{x, C + hi>)\ > > 0, j = 2, . . . ,5 for < h < h{(, v) small enough. 
Then there exists an exact solution 9{x, C, i^, h) such that for any 5^ < 13 

(2.9) 6 - eUo^U^.f-^.^)'^^ [Tei + 0{h)] < /ie~'^*^>H(?' on [a,oo), 
where T = T{x, Q is given by (I1.13p . ei = (1, 0, 0, 0, 0), and 

(2.10) = ^i(x,C + H + 0(/ie-'^^'). 

Separation of eigenvalues. Before giving the proof of Theorem 12. II we check the hypothesis 
(j2.8p for all C. Since < -R and we take h small, it suffices to check (j2.8p when v = 0. We have 



(2.11) 



Since s{x, C) = \/ C'^ + CqT] and 

(2.12) K^C' + cgr? > KC, 

it follows in view of Assumption 11.31 that when 5iC > 0, hypothesis (j2.8p is satisfied on [0, oo). 
When = and |C| > sup^, corj^^"^, we have 



(2.13) \C\>\s{x,C)\ = ^J\C\^-clv>Cc>0, 

so dSS]) holds on [0,oo). When = and |C| < inf^ cor/^/^ then 



(2.14) |Cl>s = Vco^-|CP>Cc>0, 

so again (j2.8p holds on [0, oo). 

Finally, consider ( G Class III \ {±icor?^''^(oo)}. In cases M and D ()2.13p holds on [x^, + 5, oo) 
(where x^, denotes the rightmost turning point when there are two turning points corresponding to 
^ in case M), while in case I (12.14p holds on [x* + (5, oo). Thus, for ^ in Class III (12. 8p holds on 
[x^, + (5, oo). 

Remark 2.2. 1) Observe that for in a compact subset K C {{^C ^ 0} \ Class III) the constants 
in ()2.13p . ()2.14p can be replaced by a uniform constant Ck > 0. Thus, for C € -ftT and < R, 
the separation inequality ()2.8p holds with uniform constants C^- and h{K,R) on [0,oo). 

Similarly, one obtains uniform constants Ck, h(K,R) in (j2.8p on [x.^ + 5, oo) for in a 
compact set K of the allowed set of Class III values for each type of profile. 
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Proof of Theorem \2.1[ 1. We change to x coordinates {x = xh, h = i), using (•) to denote ^ and 
(') to denote The spectral separation (I2.8p allows us to define a smooth conjugator Si{xh, C, f, h) 
such that 

(2.15) S^'i<l>oiS:hX) + hMxh,i^))Si = (^^^ •= ^i' 

where /i^ is given in ()2.6p . and the eigenvalues of Gi are j = 2, . . . , 5 as in (12. 6p . We can take 
the first column of Si to satisfy 

(2.16) Siei = Tei + 0{h). 
Since the eigenvalues of Gi satisfy 

(2.17) JR/i* > 0{he-^''^) 

we can find another smooth conjugator S2{xh,(^,h',h) = (\ ^] such that 

(2.18) 52-1^1^2 =(f Gs)^^^, 
where 

(2.19) KG2 > 0(/ie-'^''*). 

Set 53 = SiS2- Then the function W defined hy 9 = S3W satisfies 

(2.20) W = {02- hS^^S!i)W. 

Applying (12. Sp again, we can find a conjugator of the form Si^xh, C, z^, /i) = / + hS/^ such that Y 
defined hy W = S4Y satisfies 



(2.21) Y 
where B = 0{e~^^^) and 



"0 cj-^'^^- 



3f?G3 = RG2 + Oihe-l^^h > 0(/ie-^^*) 

(2 22) 

^5 =/it + 0(/ie-^*^). 



2. Next we set Y = g/o" A'U^^.Ci',^)^^^- ^nd observe that Z satisfies 



(2.23) Z={^ ^ ^^)Z-h^BZ, 





Gs-fi\I 

where 

(2.24) K(G3 - > 0{he-'^'^^') 
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We can now apply Lemma lA. II (with any number less than /3) to obtain a solution Z satisfying 



(2.25) |Z-ei| < C/ie-'^**'' 
by ()A.5p . Tracing back we find 

(2.26) \Y - eio''^i(^'''---)'^^ei| < C'/ie"'^*^'''|e'^o'^?(^'^'-)'^^|, 
so 

(2.27) \W - elo [ei + hS4ei] \ < C^he-^*^^\e-fo ^''^^'^'-^'^'l. 
Since 9 = S3W and S^ei = Tei, 

(2.28) \9 - e/o^/^Us/^vOrfs ^ q^;^)] | < (j^j^^-S^xh^J^ f,l(sh,...)d-s^ [^,00). 



Switching back to x coordinates gives (with s = sh) 

(2.29) \e -eT^^^ ^'^li'^■■■)'i'[Tel+0{h)] \ < C^he-^*''\eTi ^^^^'^■■■^'^'\ on [a, 00). 

□ 

As an immediate corollary of Theorem 12.11 we have 
Corollary 2.3. When C, ^ Class III, the exact decaying solution 9{x,T,e) of (jl.7p satisfies 

(2.30) \0{0,T,e) - r(0,C)ei| < C^/e for e > e(C) sufficiently large , 

where r = + i'. Thus, the first column of the matrix (|1.13p . evaluated at x = 0, gives the desired 
approximation. 

For Q in a compact subset K C ({KC > 0} \ Class III) and < R the constants C\, e{() in 
()2.30p can be replaced by uniform constants Ck, e{K,R). 

Proof. For such the formula of Theorem 12.11 applies on [0, 00). The uniformity statement follows 
from part (1) of Remark 12. 2^ the uniform choice of ho in Lemma lA.H and inspection of the proof 
of Theorem 12.11 □ 

Remark 2.4. Theorem 12.11 and its corollary do not require real-analyticity of the profile w, and 
are true without the restriction that w be of type I, D, or M. 

2.2 Comparison of exact and approximate solutions for ( E Class III. 

When C lies in Class III, a matching argument is needed to analytically continue the exact 
solution 9 on [x* + 6, 00), or an appropriate multiple thereof, around the turning point x*. To carry 
this out we must first examine how 9 compares to the approximate solution 9i = Tvri ()1.15p for 

X > x^, + 6. 

The expression for 9i is given in (jl.l5p - (|1.17p 

(2.31) 9i{x,T,e) = e^/o'^i(^''^)'^^+i'o:c„^ii(«'C'-)'^- [T{x,C)ei + 0{h)] . 

We assume C £ Class In case D (resp. I) we take Cw = C_ (resp. C+), a contour on the 

real axis except for a short excursion around in the lower (resp. upper) half plane. In the En 
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integral for case M, a turning point where d > is excised via the upper half plane and one where 
fi < via the lower half planei^^l 

Cases D and M. In case D the formula (|2.29p for 9 holds on [x^ + 5, oo) where 

(2.32) n{ = mix, C + uh) + 0{he-'^'=), /3 > 0. 

With fii = fii(xX) and fi\ as in ()2.32p here and in the rest of this section, we may write 

(2.33) e'^-'fof^he, = e""^'>^+^oV E,.+h-^ f,^(^{~,,)-f,^ E,,^^^ ^ e^^(^'^''^)0i(l + 0{h)) 
where 

(2.34) H,{x, C, u) := e /"(/i" - " T En- 

Jo Jo,C- 

Thus, ([2:29]1 and ([233]) imply 

(2.35) 1^ - e^'^^'^'^^^^il < ^|e^^(^''^''^)^i| for x > x, + S. 

The inequality (j2.35p holds also in case M, where x* now denotes the rightmost turning point 
in cases where there are two turning points. We do not treat the special value C = icori^^'^{xM), 
where xm is the location of the maximum. 

Next we examine the growth rate of e^^ . 

Proposition 2.5. For profiles of type D and C, G Class III-\- \ {cory^/^(oo)}, the factor e^^(^'^'^) 
satisfies the estimate: 

(2.36) Cie^('^i'^)^ < |e-^^(^''^''^)| < Cse^^"^^)^ for x £ [0,oo), x + x„ 

for some positive constants Ci, ai depending on and the profile w, but which can be taken 
independent of e, and v, and a compact subset of C, as above. For "^v > and > (resp. 
9u < Oj (j2.36p represents exponential decay (resp. growth) in x. 

For profiles of type M, we also exclude the special value (" = icor]^^'^{xM)- Then the estimate 
(j2.36p holds away from turning points with uniform constants for a compact set of C,. 



Proof. 1. Since we have 

-K ('k(C + hv) + Y^(C + hvY + Co?? 
5i/ii(x, C + hv) = ^ ^ ^ < 

(2.37) 



mu — mc + hu) = -{ — — mu. 

r]u r]u \ rju J 

For X < x^: -\- 6, X ^ x^,, it follows from (j2.37p and the formula for En ()A.24p that there exist 
positive constants Ci as described such that 

(2.38) < le-^^^^'^''^)] < Ca. 



^^When there are two turning points X2t < xi, in case M, starting with s = i\s\ in a; < X2*, analytic continuation 
along such a contour gives s = \s\ on X2, < x < xit and s = i\s\ on x > xi,. 
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2. Now assume x > x^, + 6. For such x we write 

H^{x, C, u) = H^{x:^ + 5, C, 1^) + H^^{x, C> z^) where 



1 /"^^ 



(2-39) ..b. . ^ 1 

Using (j2.37p and C) = (since t(; is of type D), we see that 

(2.40) 3? - fii) - Eu) = -—^u - — iKJ(C + H^ + cgf/ - K^ii + ©(e" 

\n J rju rju n ^ 

Assumptions HTBt 11.41 and the formula for En ()A.24p imply (with s = y^C^ + CqT? 



2 A 

(2.41) En{x, C, = -—y - —-y + 0(e-^"). 

•qu rju s 

The function s is bounded away from for x > x* + (5, so we can expand the square root in (j2.40p 
about ('^ + Cgr/ = obtaining: 

(2.42) -^(C + huf + cl^i = -s{l + + ^ - ...). 
Since s = i\s\, we get after noting cancellations 

(2.43) ^ - fii) - En^ = ^-^5R(z/iz.2) + 0{h^) + 0{e-P-) 

for X > Xtf + 5. In view of (|2.38p . (|2.39p . and Assumption 11.31 this implies the estimate (j2.36p with 
possibly new Cj. 

3. With C = icQVi^l'^{xM) excluded and redefined as described before the Theorem, the 
proof for case M is the same as above. 

□ 

Remark 2.6. For rate functions with r^, 7^ at equilibrium and satisfying the weaker assumption 
discussed in Remark II. 5 ^ we see from the formula for En that there will be an extra non-decaying 
term in (j2.4ip depending on r^. However, in ()2.43p that term will be cancelled by a corresponding 
extra term in the expansion of yuj. 

Case I. In this case we have s(x, C) = > C > for x > x* + 5, and computations almost 
identical to those above show that 

(2.44) 3? - - En^ = ^J^{hu^) + 0{h'') + 0(6"^^) 
for X > x^K + 5. 

Proposition 2.7. For profiles of type I and C, G Class IlT\{cQr]^^'^ {00)} , the factor e^'^^''^''"^ satisfies 
the estimate: 

(2.45) Cie^i-^^4> < |e^<(^'C''^)| < C2e'^(-»^4> for x e [0,oo), x / x,, 

for some positive constants Ci, ai depending on x* and the profile w, hut which can he taken 
independent of (j^, e) and a compact suhset of C ahove. For Kz^ > and < Jftzv (resp. 
> '^v > (I2.45P represents exponential decay (resp. growth) in x. 

Remark 2.8. In cases I, D, and M the estimate (|2.35p implies that the statement M(€,C,v)6 ~ 61 
("~" as defined below (11.24p ) is not true for any multiple M(e,C,z^) of 9 on [x* + 5,b) for any 
6 G (x * +5, 00]. 
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3 Exact solutions on bounded regions near a turning point 



This section is devoted to proving Theorem 13 -H which provides a sufficient condition for approx- 
imate solutions 9i of 9' = —P*6, defined by formulas like (|1.15|) . to be close to exact solutions when 
e is large. Theorem 13. II was informally stated in [Er3| with a reference to a set of 1954 N.Y.U. Notes 
by K. O. Friedrichs. The result is referred to in [Er3] as the "method of the parameter problem". 
The proof we give here is an adaptation to ODEs on bounded regions of the complex plane of a 
similar argument in Chapter 6 of |CLj for ODEs on bounded intervals of the real line. 

3.1 Method of the parameter problem. 

Consider a general N x N system with parameter e >> 10 

(3.1) 9' = e^{z,e)e 

on some bounded, simply connected region Ti C C For i = 1, . . . , N and m > fixed, let 

(3.2) Oi = e^^»(^)+^'»(-) [bio + + ■■■ + e-^^^+^^i^m+i)] := e'^'^'^'^ Pi{z, e) 
and define N x N matrices 

(3.3) Q = dmg{qi,...,qN), P = {Pi, . . . , Pn}- 

We assume that P~^{z, e) is uniformly bounded for z G Ti, e large and that the 9i are approximate 
solutions of (13.11) in H of order m in the sense tha10 



(3.4) (Pe^)' = (P' + PQ')e^ — e^mi^e^ = e^Pe^ + e($^ - ^)Pe^, 
where 

(3.5) \^-^m\<^ior z£n, 

and thus also |e($ - ^rn)Pe'^\ < C-j^ for z £7i. 

The following Theorem gives conditions under which we can find an exact solution 9 of (|3.ip 
close to 9i in the sense that for e large and z G Ti: 

(3.6) \0-0^\<cJ—^mn. 

Of course the labeling of 9j is arbitrary, hence the result applies equally to 9j, j ^ 1. 



Theorem 3.1. Consider the N x N system (j3.ip with parameter e >> 1 on the bounded, simply 
connected region H C C Let %' be a bounded, simply connected neighborhood of%, and suppose 
that <I> is analytic for z £ T-L' . Suppose that the approximate solutions 9i given by ()3.2p are analytic 
in Ti' and satisfy ()3.5p for some m> 0, and that P~^{z,e) is uniformly bounded for z £ Ti' and e 
large. 



^^Here we are for convenience suppressing dependence of $ on additional parameters in the notation. All the 
arguments below work when there is continuous dependence on a compact set of parameters. 



i^Thus, = HP' + PQ')P-^. 
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a.) For each j G {2, . . . , A^} suppose there exists a point Zj E T-L' such that any point z €z Ti can 
be joined to Zj by a path lying in %' on which ^{hj — /ii)(t, e) decreases (not necessarily strictly) as 
T moves from Zj towards z. For any z £ H let Vj{z) be such a path starting from Zj and suppose 



N 



(3.7) V sup \Vj{z)\ <L<oo, 

where \'Pj{z)\ denotes the length ofVj{z). Then there exists an exact solution 9 of in % such 
that (|3.6p holds. 

b.) There is an obvious analogue of part (a) (not requiring analyticity) for the case where % 
is taken to be a bounded closed interval [a, b] on the real axis and the paths Vj are chosen inside 
that interval. In this case existence of admissible paths is equivalent to the assumption that for any 
fixed j 7^ 1, either ?fifj,j < K/ii or else Rfij > on all of [a, b], i.e., there is a neutral spectral gap 
between fii and all other jij. 

Proof. For convenience we take = 3 in the proof. Pick Z2 and z^ satisfying the assumptions in 
Theorem 13.11 and define 3x3 matrices 

(3.8) P^^'^ = [Pi,P2M. P(') = [0,0,P3] 
satisfying 

(3.9) p = pm^p(?.)_ 

Observe that 9i is a solution of 9' = e(^„fi and that (|3.ip may be written as 

(3.10) 9' = e^mO + e{^-^m)0. 

We construct a solution of (j3.10p in as a solution 9{z,e) of the integral equation 

9{z,e) = 0i(z,e) + eP(2)(^,e)e<3(^'^) f e-Q(-.^) [p-i(r,e) (<I>(r,e) - f„^(r,e))] 9{T,e)dT+ 

(3.11) ^ 

eP(3) e)e«(^'^) r e-«(^'^) [p-i(r, e) ($(r, e) - $™(r, e))] 9{t, e)dr, 

where the first and second integrals are on paths V2{z), V^{z), respectively. Note that the only 
exponentials appearing in the elements of e)e'^*^^'^-'~'^*^'^''^) are e*^^''^^~*^'^'^\ j = 1,2, and the 

only exponential appearing in the elements of p(3)(2:, e)e'5'^^'^)~'5('^''^) is e''^^^''^^"'?^*^'^''^^. 
We solve ()3.1ip by the iteration scheme 

(3.12) 

0(,+i)(z,e) = 0i(z,6) + 6P(2)(^,g)eQ{^,^) j\~Q{r,e) [p-\r,e) (<I>(t,6) - 1>™(t,6))] e(o(r,6)dr+ 
where 0(o) = 0- We claim 

(3.13) KVi) - < ^ sup |(0(,) - 0(,_i))e-'?^(-'^)|, 
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m > 0, for C2 depending on the norms of P, P~^, and Ci (from (|3.5p ). as well as the Kfej. In 
estimating the contribution to the right side of (j3.13p coming from the first integral in (j3.12p . for 
example, we have used (j3.5p . (|3.7p . the uniform boundedness of P~^, and the fact that 

(3.14) \gg2{z,e)-g2{T,e)-giiz,e)+gi{T,e)^ _ ^n{g2iz,e)-g2(T,e)-gi{z,e)+gi{T,e)) ^ ^ ^ 'P2{z), 

where M2 depends only on "Siki, i = 1,2. Thus, for e large we obtain a contraction and the 
converge uniformly on ^ to a solution of the integral equation. Moreover, if |^ie~''^| < Co, then 
for e large enough p.l3p implies 

(3.15) 1^(^)6-^1 (^'^) I <2Co. 
From the integral equation we then obtain 

(3.16) \{6-9^)e~'^^(^'^^\<C2 

□ 



Definition 3.2. a) Let 6i he an approximate solution of ()3.ip satisfying (j3.4p for some m > in 
a region Ti C C We say that 9i is an approximate solution of order m in H. 

b) Suppose that 9i as above has been shown to be close to an exact solution on some subregion 
Ha C H in the sense of ()3.6p . Then we say that 9i is admissible in Ha- The curves defining the 
boundary of the admissible subregion Ha o^f^ often referred to as Stokes curves. 

We will see below that the admissible region T-La is sometimes a proper subregion of H. For 
example, in applying Theorem 13.11 to justify an approximate solution 9i in the case where A'^ = 3, 
it may happen that one can choose paths V2{z) satisfying the requirements of the Theorem for all 
z GTi, but that one can choose paths 7^3(2;) only for z in a proper subregion Tia of H. In that case 
9i is admissible only in Tia- 

Definition 3.3. When applying Theorem \3.1\ to test the admissibility of 9i, we refer to paths 
on which ?R.{hj — hi) decreases as Vj paths for j G {2,...,A^}. More generally, when testing 
the admissibility of 9k, we refer to paths on which ?fi{hj — h^) decreases as Vj paths for j G 
{1, . . . , A^} \ {k}. The context (i.e., choice of 9^ being tested) should prevent confusion. 

Remark 3.4. Consider again the 5x5 system 

(3.17) 9' = -P\x, T, e) = [e^o{x, Q + <i>i(x, u)]9 

with /Lij, j = 1, . . . , 5 as in p.lOp . When applying Theorem 13.11 to test the admissibility of 6*1, for 
example, observe that since ^3 = = //s, the paths, P4 paths, and "Ps paths can all be taken 
to be the same. In this situation we usually refer to all of these paths as paths. 



3.2 Analytic continuation and the choice of integration paths. 

We now fix (" = i\C,\ G Class III. By Assumptions [L2] and [L3l the coefficients in the system 
(j3.17p extend analytically to a neighborhood of [0, 00). Quantities appearing in the definition of 9i 
()1.15p are of two types: those that depend real-analytically on the profile and those like ^i, ^2, 
and certain terms in T{x, () (jl.l3p and En ()A.24p that depend on s{x, (^) = y^C^ + Cgr/(x). Terms 
of the first type extend analytically to a neighborhood of [0, 00), while analytic extensions of terms 
of the second type have branch points at turning points where s(x*, () = 0. 
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When testing the admissibility of a solution 9i in a region T-L whose points all lie close to a 
turning point where CqT] is either increasing or decreasing, there is a simple procedure given 
in |Er3j for choosing paths that satisfy the requirements of Theorem 13.11 We need to use this 
procedure in section HI so we illustrate it here for the reader's convenience in the increasing case 
where 

(3.18) d := ^^(x,) > 0. 

ax 

We take Ti to be an open disk centered at from which a small (closed) neighborhood of 
the segment [0, x^] has been removed, and explain first how to choose V2 paths and V3 paths for 
showing the admissibility of 9i. Suppose that 9i is defined before continuation into T-L by formula 
(|23T]1 for X > X,. 

Expanding s'^{z) = + (^oVi^) about x^,, we compute from the formulas (|2.1ip : 



(3.19) 



H2iz) - fii{z) = {2Kd'2 /r]u){x^){z - x*) 2 [1 + 0{z 
Hsiz) - fii{z) = (i/??^K)(x*)[l + 0{{z - x*)5], 
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where {z — x*)2 denotes the branch that is positive for z > x^. We set 

(3.20) hji = hj - hi, h*i = hji{x^), 

and integrate (j3.19p to obtain 

h2i{z) - h*2^ = {4Kd^2/3r,u){x^){z - x,)t + 0{{z - x)t) 
(3.zlj jj 
h3i{z) - = {i/r]2K){x^){z - X*) + 0{{z - x*)2). 

Before choosing paths on which K(/iji — h*^) (and therefore also K/iji) is nonincreasing, we first 
locate the transition rays where K(/iji — h*-^^) ~ 0, j = 2,3. With (p = arg(z — x*) these rays are 

(3.22) (P2i=±^,7r; 03i=O,7r, 

where the subscripts have the obvious meaning. The transition rays defined by (p2i determine 3 
sectors centered at x*: 

(3.23) So:-^<(t><^, Si:^<(t><TT, 5_i : -vr < -|. 

Paths in T-L on which 3f?(/i2i — /i2i) nonincreasing are easy to determine if one knows the level 
curves of 3f?(/i2i — ^21)- sector Sq, 5i(/i2i — /i2i) is positive with level curves that are roughly 
parallel to the union of the two rays </> = f and = — ^. In sectors S±i K(/i2i — /i2i) is negative 
and the level curves have a similar description. V2 paths with initial point Z2 > x* can now be 
constructed by piecing together segments that lie on level curves of 3?(/i2i — /121) with (correctly 
oriented) segments that are transverse to level curves. If T-L is small enough, any point z & Ti can 
be reached by such a path starting from Z2- 

From (j3.2ip and (j3.22p we see that the level curves of 3f?(/i3i — h'^^) are well-approximated by 
horizontal lines near x* and K(/i3i — h^i) > in Q{z — x*) < 0. Thus, if we choose Z3 near x* in 
$5(2: — X*) < and if T-L is small enough, any point z G Ti can be reached by a path starting at 

Z3- 
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A similar procedure applies when testing the admissibility of approximate solutions j = 
2, ... ,5. For example, when testing 62 in T-L one chooses Vi paths and V3 paths usinel^ 

1 3 5 

huiz) - = -{iKd^ /3r]u){x^){z - x^)2 + 0{{z - x^)2) 
(3.24) ^ 3 

h32{z) - /l32 = {i/r]2K){x^){z - x^) + 0{{z - X*)2). 

4 Matching arguments for ( G Class III 

In Theorem 12.11 we constructed for ( € Class III an exact decaying solution 9 of equation (12. 7p 
on [Xif +5,00). In this section we give the matching arguments needed to determine the analytic 
continuation of that solution up to a; = 0. 

4.1 Case I: Cq — is strictly increasing. 

In order to carry out the matching argument we must first use Theorem 13.11 to construct a basis 
of exact solutions on bounded domains near x^:. 

Let ^j(x, C, i^, e), i = 1,...,5, be approximate solutions constructed in the standard way and 
defined before analytic continuation for a; > x=k by formulas 

(4.1) 0i{x, C, V, e) = exp [ehi{x, C) + h{x, C, i')]ai{x, C, 1^), 

where 

/•X 

hi{x,()= / /ii(r, C)dT, for alH, 



(4.2) 

ki{x,C,^) = / Eii{T,C,u)dT, fori = 1,2, 
Jo,c+ 

and the subscript C4. indicates a contour on the real axis except for a small excursion in the upper 
half plane to avoid x=k In (j4.2p the square root s(x,C) = i\s\ for x < x*, s = \s\ for x > x^, and 
is given by analytic continuation on C+. The functions En and E22 have terms with (z — x*)"-' 
singularities at x*, where j G {1/2, 1}. On the other hand /ij is bounded near x^. 

Definition 4.1. For 6 > small and any set S C C, we define a 5— neighborhood of S to be the 
union of open balls of radius 6 centered at points of S. 

Below we let denote the closure of a 5— neighborhood of the ray {z : arg(z — x^,) = (p}. Let 
B C C denote a 45-neighborhood of the segment [x* — 66. x^]l^ 

Let 9i^M be an approximate solution defined by analytically continuing 9i to B\Nt^, and denote 
by 9i^M the corresponding exact solution justified by Theorem 13.11 in this region. Paths of type V2 
in the sense of Theorem 13.11 can be chosen starting at some Z2 > x^, while V3 paths can be chosen 
starting just below x=k in < 0. Here one can use exactly the same paths as those constructed in 
section [3^2] using (I3.2ip . 

Denote by ^2 an exact solution constructed by a similar procedure, but based on the approximate 
solution 92 analytically continued to B\ N^/^. Paths of type Vi can be chosen starting at a point 
Zi in the 3rd quadrant (relative to x*), while V3 paths can be chosen starting in 9z < 0. 



^^The choice of branch of {z — a;*) 2 in these formulas depends both on the region H and on the specification of 62 
before its analytic continuation into H. Typically, 9j is specified first on either x < x, ot x > x*. 

^''For the matching arguments to be valid we need to work near x,, so the parameter 5 introduced here is to be 
fixed sufSciently small. For the continuation of solutions all the way to a; = see step 4 of the proof of Proposition 
K3\ 
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Similarly, we let 6*3, ^4, ^5 be exact solutions associated to analytic continuations of 03,64, 9^ and 
defined in i? \ A^_^/3. In this case paths of types Vi and can be chosen to start directly above 
X* on the axis ?R.z = x^. 

Let 9f denote an approximate solution obtained by analytically continuing 9i to B \ A^_7r/3 
and let 9^ denote an an approximate solution obtained by analytically continuing Oi to B\ N^/^. 
Each of these solutions can be justified by Theorem 13.11 on a 2(5-neighborhood O2S of the segment 
B n [0, X* — 36] , yielding corresponding exact solutions 9f and 9^ . For Of the paths of type V2 and 
V3 start in the third quadrant, while for 9^ paths of type V2 start in the second quadrant and 
paths start in the third quadrant. 

It is helpful for computing F(r, e) to have explicit formulas for 9f{z,C,i',e), z G O^?,- These 
have the form 

(4.3) 9f{z, C, V, e) = exp \ehf{z, C) + kf{z, (, iy)]af{z, (, u) 
where 

(4.4) ehf{z,C) + kf{z,C,u):= r [efii{T,C) + Eu{T,C,u)]dT 

Jo,p± 

and af are obtained from T{x,C)ei, x > x*, by analytic continuation. The path V+ consists of 
the segment [0, — 26] followed by any path in O25 from x^, — 26 to z, and the function of z given 
by s{z,() = \/ C"^ + <^Qr]{z) appearing in fii is obtained from the branch taking positive values on 
X > x^, by analytic continuation in B^ \ N_^/3. The same square root, call it s+{z,() is used, of 
course, wherever it appears in af. It is readily checked that in O2S the function 9f defined by (|4.3p 
equals 9f as originally defined. 

The path P_ = V^^a^'P~,b, where V-,a is a segment from to x* — 26, followed by a semicircle 
in Qz > terminating at x* + 26. The path V-^b starts at x* + 26, continues on the semicircle in 
Qz < terminating at x* — 26, and continues from x* — 26 along any path in O2S ending at z. In the 
part of the integral (14. 4p on V^^a the square root s+(z, () is used, while on V-^b, one uses S-{z, (), 
which is obtained from the branch taking positive values on x > x,, by analytic continuation in 
B+ \ N^/^. The square root s_(z,C) is used wherever s appears in . In O2S the function 9^ 
defined by (I4.3P equals 9^ as originally defined. 

Finally, observe that since //3(x,C) = ^, the function /i3(x,C) extends analytically to a neigh- 
borhood of [0,cxd) containing B. 

Remark 4.2. The approximate solutions 9f defined here in i? \ -^=p| extend analytically (as 
approximate solutions) to a neighborhood of the segment [0, x^, — 26], and explicit formulas for 
those extensions are given by the formulas (I4.3|) . (l4.4p using the obvious prolongations of the paths 
V±. The admissibility of the extensions on [0, x* — 26] is shown in step 4 of the proof of Proposition 

Proposition 4.3. For a profile of type I, fixC, ^ Class ///+\{«co??^/^(0-|-), ico?/^^^(oo)}. Let 9 he the 
exact decaying solution on x > x.^ + 6 constructed for case I in Theorem \2.1\ extended analytically 
to [0,00). For an appropriately selected multiple 9i^g = G{e,C,,u)9 we have on [0, x* — 2>6] 

(4.5) 9i,G = 0t + 9^ + 0{l/e). 

Here 9f and 9^ are analytic extensions of the exact solutions 9f defined above on O25 which satisfy 

(4.6) ]9f - 9f] < -^]9f] on [0, x* - 3(5], 
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where 9^ are given by the formulas (j4.3p . (j4.4p . 

The constant in (|4.6p can be replaced by Ck for C, in a compact subset K of the allowed 
values. 



Proof. 1. Let xi = x^, -\- 25 and recall (j2.35p 

^ I , n E Ml + (n -El . 

ai. 



(4.7) |0_e^^(-'C''^)0i| < %H.i-^^^u)Q^\ ^^^^^ ^ g^/o>i+/o!c+^ii 

e 

Multiply through by 

(4.8) G(e,C,z^) :=e-^^(^i'^''^) 
to get 

(4.9) |G(e,C,i^)^-e-^^(^'^'^)-^^(^i'^'^)6'i| < -\e"^i^^'^^'^)-"^i^^Meil 

e 

Evaluating this at xi we obtain 

(4.10) \G{eX,^)0-ei\{xi) < ^\ei\{xi), 
2. Next we carry out the matching argument using 

(4.11) ei,G{x,C,i^,e):=G{e,C,iy)e 

for the exact decaying solution on [xi, oo). Observe that the exact solutions 6i^g and 9i^m, ^2,- • • , ^5, 
have original domains (i.e., before analytic continuation) that include [x* + 2(5, x* +4(5)o Using 
the basis {Oi^m, 02, ■ ■ ■ , db} we have 

(4.12) 9i,G = CiOiM +0202 + ■■■ + C5^5, 

where for all i, Ci = Ci(e, i^)- By Theorem 13.11 we have 

C 

(4.13) \9i^M - ^i| < -l^il on [x* + 2(5, X* +45) 

and analogous statements hold for the other 0i. An estimate like (j4.13p is not true for Oi^G: but at 
xi we have from (j4.10p 

(4.14) |^i,G-^i|(xi) < ^l^il(xi). 

The coefficients in (I4.12p are independent of x so we can evaluate equation (|4.12p at xi to 
determine them. Recall Rhs^xi) = and let 

(4.15) a := -K/ii(xi) > 0, 6 = K/i2(xi) > 0. 

A direct analysi^^ of the coefficients in ()4.12p using ()4.13p and (j4.14p shows that 

(4.16) ci = l + 0(e-i), c2 = 0(e-(''+^)7e), cs = ©(e-'^Ve), cs = 0(e-'^7e). 



"^^AU of these exact solutions can be analytically continued as solutions to a neighborhood of [0, oo), but it is only 
on the original domains that we have information about asymptotic behavior as e — ^ oo. The original domain of 
9i,M, for example, is _B \ A''^. 

^^This statement about the coefBcients is proved by writing 9i,g = e'''^^'^''^ai + 0{e'''^^^'^^ /e) and similar expressions 
for the other solutions, plugging into equation (|4.12|) . using the linear independence of ai, . . . , as, and taking account 
of the growth rates of the exponentials. The straightforward details are omitted. 
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The functions of x in (j4.12p extend analytically to a neighborhood of [0, oo), and so (j4.12p holds 
on that neighborhood. 

3. From the earlier construction we already know the asymptotic behavior as e — )• oo of 
02, ■ ■ ■ ,0^ for B (1 {x < x^}, but that is not the case for 9i^m, whose behavior we know only on 
B\Nj^. To remedy this we expand Oi^m in terms of the basis {Of , 9^ ,9^, . . . , 9^} on the intersection 
of the original domains: 

(4.17) ^■=(.B\ N^) n O25 n {B \ N_^/3). 
This gives 

(4.18) 9i,M = c+9t + c^9^ + ^3^3 + • • • + d^h, 

where again all the coefficients depend on (e, C,i^)- Next we evaluate equation (j4.18p at a point 
z £lr\{Qz < 0}, and use ([O]), ([MI]), and h^^ = to see that 

(4.19) a+ := K/i^(i) < < := K/i^(l) < 03 := ^h^iz). 
We have 

(4.20) Cie'^t' < \9t{z,C,iy,e)\ < Cse"^, 

and estimates analogous to (|4.20p at z_ hold for {^i,A/,^f} and {^3,^4,^5} with af replaced by 
and 03 respectively. 
Using 

(4.21) \0^^^^-O^\(^)<^\e^\{z) 

and an analysis of coefficients similar to the one that gave (|4.16p . we find 

(4.22) c_ = l + 0(l/e), |c+| < Ce^"^-"^^', dj = 0(6^"^ = 3,4,5. 

To get more information on c+ we select G Xn {9z > 0}, evaluate equation (|4.18p at z^, and 
define 

(4.23) := m];{z^) < < bf := mf{z^), and 63 := ^hiz^) < 0. 

Now {9i^M,9t} grow like e^i^ at z'^ (in the sense of (|4.20p ). 9^ grows like e^i^% and {^3,^4,^5} 
grow like e^^''. Thus, when equation (|4.18p is evaluated at z", all terms to the right of C-^-9f decay 
exponentially fast in e. So we conclude from 

(4.24) |^,^^,_^+|(^tt)<^|0+|(^tt) 

that c+ = 1 + 0(i). Substituting ()i38]) into KW we obtain 

(4.25) 9i^G = (1 + 0(l/e)) 9+ + (1 + 0(l/e)) 9^ + + ■ ■ ■ + hh, 

where the kj decay exponentially fast in e. 

4. At this point we know the asymptotic behavior of the functions on the right in (j4.25p only 
near x*. For example we know 

(4.26) l^f -6'fl < ^l^fl 
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only in O2S, a 25— neighborhood of the segment B n [0, — 35] = (x^, — 105, x^, — 35]. To conclude 
that (|4.26p holds on all of [0, — 35] we observe that Of can be constructed by Theorem 13.11 on 
[0, x^, — 55] (say) using V2 and paths approaching zero along the x— axis, but which start at the 
same initial points already chosen for the justification on 025- This works because ?R.^j = on 
[0,x=K — 55] for j = 1,2,3. The other functions on the right in (|4.25|) are treated similarly. 

5. The ability to choose (7^ uniformly for C in a compact subset of the allowed Class III values 
follows from the corresponding uniformity statements in part (2) of Remark l2.2l and the comparison 
proposition Prop. 12. 7^ and inspection of the proof of Theorem 13.11 together with steps 1-4 above. 

□ 

A simpler formula for 9^ and 9^ on [0,a;* — 35]. It follows directly from the definition of 
9^ given earlier that 



(4.27) 



K C, e) = exp 



Jo 



where s(x,C) = ^l^l whenever it occurs in (j4.27p . Let us define 



(4.28) 



9^{x,C,J^,e) = exp 



/ (e^2(r, 0+ ^22 (r,C,i^))dT 
Jo 



ai{x,(,u) on [0, X* — 35] , 



a2(x, C, J^) on [0, x^, — 35] 



using the same s{x, (); here aj = T(x, C)^ji j = 1; 2. Observe that by changing the sign of s in the 
formula for fii we obtain ^2', similarly, we obtain a2 from oi and -E22 from En by changing the 
sign of s (see section IA.4P . It follows immediately from this and the formulas for 9^ given in (j4.3p , 
(OD that 

9^ = a(e, C, 1^)02 [0) X* ~ 35]. 



(4.29) 9-^ 
Here the factor a is given by 

(4.30) a(e, C, z/) := exp 



/ (e^i + Sii)dr 
Jc 



where C proceeds from along the x— axis to — 25, then around x* on a clockwise circle of radius 
25, then back to along the x— axis. The square root changes by analytic continuation along the 
contour, starting on [0,x* — 25] with s = i\s\ as abovj^. 
Defining the exact solution 92 '■= a~^9Y , we have: 

Corollary 4.4. For a profile of type I, fixC £ Class ///+\{icor/-'^/^(0+), icor/^/^(oo)}. On [0, x*— 35] 
(4.31) 9i,G{x,C,^,e) = 9+ + a{e,C,iy)9+ + 0{l/e), 

where a is given by (|4.30p and 



< 



j = 1,2 



(4.32) 

on [0,x* — 35]. In particular, we have 

(4.33) |^i,g(0, C, e) - [r(0, C)ei + a(e, C, z^)T(0, C)e2] | < CJe 

for T as in ()1.13p . The constant can he replaced by Ck for C, in a compact subset K of the 
allowed values. 

Remark 4.5. Since ^j, i = 1, 2 are purely imaginary on [0, x* — 35], the functions 9~j , j = 1,2 are 
0(1) there as e — t- 00. Thus, the error term in ()4.3ip is small compared to the other terms. 



On the return trip from Xt — 2(5 to we have s 
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4.2 Case D: Cg — is strictly decreasing. 

To construct the exact solutions needed for the matching argument in this case we start with 
approximate solutions 6i given by 

(4.34) 6'i(x, C, V, e) = exp [ehi{x, Q + h{x, C, v)]ai{x, C, 

but now they are defined before analytic continuation for x < x-t where 

1 fii{T, C)dT, for all i 
Jo 

= / Eii{T,C,u)dT, fori = 1,2, 
Jo 



hi{x,C) = 

(4.35) 

ki{x,C,i' 



and s{x, C) = |s| for x < x* 1^1 

Again let i? C C denote a 45-neighborhood of the segment [x^ — 65, x,,], but now let A'^ denote 
a closed (5— neighborhood of the ray {z : arg(j;* — z) = (j)}. Let Oi^m be an approximate solution 
defined by analytically continuing to S \ N and denote by Oi^m the corresponding exact 
solution justified by Theorem 13. II in that regiono When choosing V2 paths, for example, we now 
usf@ 

(4.36) h2i{z) - h*^ = -{AK{-d)^2 /3r]u){x^){x^ - z)l + 0((x* - z)t) 

in place of (j3.24|) . Paths of type V2 and can be chosen starting in the fourth quadrant (relative 
to X*). 

Let 02 be an be an approximate solution defined by analytically continuing 02 to B \ N_^/^. 
Paths of type V3 terminating at any point in this region can be chosen with a common initial point 
in the fourth quadrant. Paths of type Vi terminating at any point in the region Bi := {z : ^ < 
arg(x^, — z) < ^} can be chosen with a common initial point in the first quadrant. Paths of type 
Vi terminating at any point in the region 82 ■= {z : — ^ < arg(2;* — z) < it} can be chosen with 
a common initial point on x < x^,. By Theorem 13.11 we obtain one exact solution 621 in Bi and 
another 622 in B2 satisfying 

(4.37) \e2j - 021 < ^1^2! in Bj, j = 1, 2. 

Similarly, we let ^3,^4,^5 be exact solutions associated to analytic continuations of ^3,^4,^5 and 
justified in B \ A'^/3. In each case paths of types Vi and can be chosen to start directly above 
x^ on the axis 3ftz = x*. 
Now the goal is to prove 

Proposition 4.6. For a profile of type D, fixQ ^ Class ///+\{ico7?^^^(0+), ico7?^^^(oo)}. Let 6 he the 
exact decaying solution on x > x* + 5 constructed for case D in Theorem \2.1{ extended analytically 
to [0,00). For an appropriately selected multiple Oi^g = G{e,C,,u)9 we have on [0, x* — 3(5] 

C 

(4.38) l^iG - ^^il < ^|6li| on[0,x*-35]. 

e 

The constant in ()4.6p can be replaced by Ck for C, in a compact subset K of the allowed 
values. 



*The value of s in x > a;*, s — i\s\, is obtained from that in a; < by continuation in SJz < 0. 
^Caution: The ray arg(a::t — z) — — tt/S, for example, lies in the second quadrant relative to 
^Formula (|4.36|) corrects a sign error in |Er3) . p. 49. 



27 



Proof. 1. As in case I we take xi = x* + 26, define 

(4.39) G(e,C,^^) =e-^^(^i'^'^), 
and note that 9i^g '■= G{e, (, v)9 satisfies 

(4.40) \h,G - ^i,m|(xi) < ^|^i,m|(xi). 

2. Using the basis {Oi^m-, ^3, ^4, ^5} we have 

(4.41) 0i^G = cihM + C2021 + ■■■ + c^h on S n {x > xi}, 
where for all i, Ci = Ci{e,C,,v). By the construction of 9i,m we have 

C 

(4.42) \ei^M - 6i,m\ < — l^i.A/l on [x* + 25, x* + 45) 

and analogous statements hold for the other functions on the right. 

3. We have an integral formula 

(4.43) Oi^m{z, C, ^, e) = exp [ehi^Miz, C) + h^uiz, C, '^)]ai,Mi^^ C, 
where 

(4.44) ehi^M{z,C)-= [ e^i(r,C)dr, 

Jo,r 

r consists of the segment [0, x* — 25] followed by any path in i? \ N-t^/s from x* — 25 to z, and the 
definitions of /ii(r, C), ki^M^z, C, v), and ai^Miz, C, v) are readily gleaned from the definition of 6i^m 
given above. We shall also need the corresponding formulas for 62 in Bi and B2 and the functions 
63,..., 9n+4 in B\ 

4. We determine the behavior of the coefficients in (j4.4ip after evaluating (j4.4ip at xi. From 
the properties of the fii on the real axis and the integral formulas for the approximate solutions, 
we easily determine that 

(4.45) a := -3f?/ii,M(2;i, C) > 0, 6 := 5ft/i2(xi, C) > 0, and K/i3(xi, C) = 0. 
Thus we have a situation like (|4.15p and the same coefficient analysis shows 

(4.46) ci = 1 + 0(e-i), C2 = ©(e'^'^+^^Ve), cj = ©(e'^^Ve), j = 3,4,5. 

5. We know the asymptotic behavior in i?n{x < x^, — 35} of all functions on the right in (j4.4ip 
except for ^21- Thus we write 

(4.47) 021 = diOi^M + (^2^22 + ^3^3 + ^4^4 + ^5^5, 

where again all coefficients depend on (e,C,v). Evaluate equation (I4.47P at a point z given by 

(4.48) arg(x=, - z) = ^ + 5, \z\ = 25, 

a point where all functions in ()4.47p are known asymptotically. The integral formulas for the 
approximate solutions yield (with paths of integration mostly restricted to the x— axis) 

(4.49) d := -K/ii,m(i, C) > 0, 6 := K/i2(i, C) > 0, and c = ^siz, () with |c| < c(5). 
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where c{6) — )• as (5 — )• 0. Direct analysis of the coefficients in (|4.47p using ()4.37p gives 

(4.50) \di\ < Ce^^+'^^'/e, da = 1 + 0{l/e), \dj\ < Ce^^-^^'/e, j > 3. 

For 6 small we have d < a,b < b, so from (|4.41|) . (|4.46!) . (|4.47|) . and (|4.5U|) we conclude 

(4.51) \0i,G - Oi,m\ < j\Oi,m\ on [0,x, - 36]. 

6. It remains to show that (14.511) holds on [0, — 35]. The problem is that, while the functions 
on the right in (j4.4ip and ()4.47p have analytic extensions to x = 0, we do not know the asymptotic 
behavior in e of those extensions. As in case I we can use Theorem 13.11 to determine this behavior 
in the case of Oi^m, because 3ft(/i2 — hi) and K(/i3 — hi) decrease as x — )• on [0,x^ — 36]. This 
allows us to choose V2 paths and V3 paths terminating at points on [0, — 35] and starting at the 
same initial points already chosen for the earlier justification of 6i^m on B D {N \ — ^)- Thus we 
obtain 

C 

(4.52) 16*1,^/ -6*1! < -j\Oi\ on [0,x* -35]. 

This method does not work for the other functions appearing in (14.4ip and (j4.47p . For example, as 
X decreases to 0, K(/ii — /12) and 3^(/i3 — /12) both increase, and therefore Theorem 13.11 cannot be 
used to show the admissibility of 622 on [0, — 35]. 

Instead we use the differential equation to estimate directly how much these functions can grow 
as X approaches in [0, x^, — 35]. We then use the estimates on the coefficients in (j4.46p and (j4.50p 
to show that in spite of this growth, (j4.5ip still holds on [0, x* — 35]. 

We illustrate this argument in the case of ^3. On [0,x* — 35] we have for some constants c,C: 

(4.53) < c < = JR^s < C < ^fi2- 

After diagonalizing the system with the transformation T(x, Q ()1.13p and setting ^3 = TW, we 
obtain with D = diag(/Ui, . . . , /U5): 

(4.54) d^\W\^ = d.^{W,W) = 2^{eDW,W) + 0{\W\^) \W\^ > e^ni\W\ - C^\W\. 
So for e large, W and therefore ^3 grows at most by a factor 

(4.55) ^^i:'-''{\^,^i\is,0+4)ds 

as X varies from x* — 35 to y, where y G [0, x* — 35]. In the expansion of 9i^g on S n [0, x=k — 36] 
the function ^3 occurs with a coefficient of 02^3 + C3. Since C3 = 0(e~"'^/e) and 

/•Xt —35 

(4.56) a > / \^fj,i\{s,C)ds 



we obtain 

(4.57) |c3^3|(y) < -e-'''e'^"*-^s^^''^'^Ods < C^e fy lR^.,{s,Ods < ^j^^ 

e e e ' 

A similar estimate for 02^3^3 shows that the contribution of ^3 on [0, x* — 35] is negligible. The 
remaining functions in the expansion of Oi^g ^-re treated similarly. Thus we obtain 

(4.58) |^i,G-^i,A/| < j\9i,m\ on [0,x, -35], 

which with ()4.52p implies the result. 

7. The final uniformity statement is proved as in step 5 of the proof of Proposition 14.31 

□ 
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Corollary 4.7. Fix ( S Class III+ \ {icori^^'^{0+),ico7]^^'^{oo)} . In case D we have 

(4.59) |^i,g(0, C, e) - r(0, C)ei| < C^/e 

for T as in (jl.lSp . The constant can be replaced by Ck for C in a compact subset K of the 
specified Class III values. 

4.3 Case M: Cq — has a single interior maximum. 

The nonvanishing oi d = -^^^^{x*) played an important role in the matching arguments for 
cases I and D. In case M denote by xm the interior point where 

(4.60) ^(xm) = 0. 

ax 

A different type of matching argument is needed for the Class III value 

(4.61) QM = icori^'\xM). 

The turning point problem for C,m is studied in [Er3j , section 6, pages 71-98, but we shall not treat 
it here. For all other non-exceptional Class III values of Q in case M, the turning point problems 
that arise can be handled using the earlier results for cases I and D, together with their proofs. 

For each such Q there are either one (x* = xi*) or two (x2* < xu) turning points . As in cases 
I and D we define Oi^g = G{€,C,i')9, where 9 is the decaying solution constructed in Theorem 12.11 
on [xu + 6, oo) and 

(4.62) G(e, C, ly) = e-^^(^i'^'^\ xi = xu + S. 

In the Ell integral in the definition of H^{x,C,,v) a turning point where d > is excised via the 
upper half plane and one where (i < via the lower half plane. 

In part (3) of Theorem 14.81 one must correctly redefine the function 9i that appears in the 
statement of Proposition 14. 6| and the functions 9^, 9f, ^J, 0^, and constant Q(e, C,z^) that appear 
in the statement of Proposition 14.31 and Corollary 14.41 

The function 9i is defined for X2* < x < xu before analytic continuation around xu by formulas 
(|4.34p . (|4.35p . where s is given by (|1.12p on the real axis, and in the ki integration path the point 
X2* is excised via the upper half plane. 

The functions 9f^ are given on X2* < x < xi* by the same function 9i before analytic continuation 
around X2* • Their continuations and the corresponding exact solutions 9f^ are then constructed as 
in case I. The functions and are then defined as in case I. The constant a(e, i^) is given by 
(I4.30p . where the contour starts at and encircles just X2*. 

Proposition 4.8. In case MletCe Class III\ {icoV^^^{0+),icor]^^^{M),icor]^/^{oo)}. 

i) Suppose there is only one turning point x* corresponding to C CL^d that -^^{x*) < 0. Then 
the asymptotic behavior of 9i^g on [0,x^, — 35] is given by Proposition \4-6\ 

ii.) Suppose there is only one turning point x^ corresponding to C and that (x^) > 0. Then 



the asymptotic behavior of 9i^g on [0,x* — 35] is given by Proposition^^ and Corollary \4.4\ 
Hi.) Suppose there are two turning points X2* < xi* with 

(4.63) ^(,,,<„<^(|rt(,,.). 

On [x2* + 3(5, xi* — 35] the asymptotic behavior of 6i^g is given by Proposition \4.6[ while on [0,x* — 
36] the asymptotic behavior of 9i^g is given by Proposition \4.3\ and Corollary \4.4\ particular, 
^i,g(0, C) e) is given again by formula (j4.33p . with a defined as explained above. 
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Proof. 1. Part (i) (resp., part (ii)) follows immediately from the fact that for x < x^:, s'^ = 
+ CqT] > (resp. < 0) just as in case D (resp. case I). 

2. In the proof of part (iii) the functions ^i,m and Oi^m are defined near xi^ as in case D, using 
the redefinition of 6i described above. The analytic extension of Oi^m to [x2* + 3(5, xu — 35] satisfies 

C 

(4.64) I ^i,A/ - ^1 1 < - 1 ^1 1 on [xa* + 35, xu - 36] . 

by the argument using Theorem 13.11 that gave ()4.52p . The redefinition of the other exact solutions 
^21, ^22, • • • needed for the matching argument near xu is similar to the redefinition of 9i^m- 
An estimate like (j4.54p is now used to show that 

C - 

(4.65) 16*1,0 - (^i,m\ ^ — I^i,m| on [x2* + 35, xi* - 36]. 

3. The redefinition of the other exact solutions 6i^m, 02, ■ ■ ■ ,05 needed for the matching argu- 
ment near X2* is similar to the redefinition of Of^ . In place of ()4.14p we can set X2 = X2* + 36 and 
use (USD, (^651) to conclude 

(4.66) |^1,G-^1|(X2) < ^|^1|(X2). 
The rest of the proof follows as in case I. 

□ 

Remark 4.9. The results of this section show that in case D (resp. case I) the exact decaying 
solution of (II. 7p decays (resp. oscillates) to the left of the turning point, and oscillates (resp. 
decays) to the right of the turning point. In case M we have oscillation to the left of X2*, followed 
by decay, followed by oscillation to the right of xi^,. In cases D and M, when ?R.iy > the oscillation 
on the unbounded interval occurs with an amplitude that decays exponentially to zero. 



5 The Instability Theorem 

In this section we present, with some simplifications and extensions, the main stability and 
instability results of |Er3] . In particular, we give the argument of |Er3) for locating unstable zeros 
of 

(5.1) V{t, e) = T6i(r, e) + iebiir, e) - 6(0, r, e) • (rht + iehy), 

where 6 is the exact decaying solution of (II. 7p on [0,oo), and the other quantities appearing in 
()5.ip are defined in section IA.3I With r = Ce + we rewrite this as 

y(r,6)=6L(6,C,i^), 

He, C, 1^) = Ch + ib2 - 0(0, C, V, e) ■ {Qht + ihy) + e'^v[bi - 0(0, C, e) • ht]. 

For IkQ > and all three types of profiles, the integrals 6i, 62 are shown in [Er3j . p. 57-64 to approach 
zero as e — )• 00. This part of the argument, which is complicated since the integrands involve 0, 
can be eliminated if one uses instead the following simpler form of the stability function derived in 
|C JLW| . Definition 4.7, and shown there to equal V: 

(5.3) y(r, e) = 0(0, C, i^, e) ■ u;(0+) - 0(0, C, e) • {rht + iehy). 
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Using (|5.3p . one finds by inspection 

(5.4) |L(e, C, i^) - Laie, C,i^)\ < Cx/e as e oo, where La := -9(0, (, ^, e) • [C,ht + ihy], 

and the constant Ck is uniform for in any compact subset 

The next proposition, which is the main step in the proof of the theorem that follows, provides 
explicit formulas for the approximate stability function L^- 

Proposition 5.1. Let ti(C) <ind t2(C) denote the first two columns of the matrix r(0,C) P-13p . 
Under Assumptions \l.Sl \L3l and \L4\ we have: 
1 )For C i Class III 

(5.5) \Laie, C, u) - Li(C)| < Cx/e, where Li(C) := -ti(C) • {(ht + ihy), 

and the jump terms ht and hy are given in section \A.^ The constant Ck is uniform for C, in any 
compact subset K C {Class IIL)^. This result holds even for profiles that are not of type I, D, or 
M. 

2) For profiles of type D and Q G Class ///+ \ {ico?7^''^(0+), ico?7^/^(oo)}, we have 

(5.6) \La{e,C,^)-Li{0\<CJe. 

3) For profiles of type I and ( G Class III^ \ {icorj^^'^ {0+) , icQrj^^'^ (oo)} , we have 

(5.7) \Laie,C,i^) - [Li{C) + a{e,C,i^)L2{C)]\ < CJe, where L-^^Q) := -t2(C) • (Cht + ihy), 

and a is given in (|4.3U|) . 

4) For profiles of type M and ( G Class IIIj^ \ {icQVi^ ^ {x m) , icQVi^ ^ {^+), icQiq^ / {oo)} , in the 
cases described by parts (i), (ii), and (Hi) of Proposition La satisfies, respectively, (15. 6p . (15. 7p . 

m- 

In parts (2) and (3) and in the three subcases of part (4), the constant can be chosen 
uniformly for Q in a compact subset of the specified set of Class III+ values. 

Proof. The result follows immediately from (j5.4p and the formulas for 6(0, v, e) given in Corollary 
12.31 Corollary 14.71 Corollary 14.41 and Proposition 14.81 part (iii). □ 

The idea now is to study V for large e by studying the functions Li(C) and Li(C)+a(e, C) 2^)-^2(C)) 
which are given explicitly in |Er3j . The functions Li ()A.19p and L2 (|A.20p are computed from the 
expressions for T(0,C)j ht, and hy already given, and a is readily computed using the expression 
for Ell (|A.24p by letting the radius of the circular part of the contour C in (14.30p shrink to zero. 
Writing = iQ, one obtains 

a{e,C,i^) = el^^''^'''\ where 

(5.8) qrj 



^^The rate of convergence to for L — La, is computed in Er3 for profiles I, D, or M to be 0{e ^^^) with no 
discussion of uniformity with respect to By (|5.3|l the rate Ck/^ holds, even for profiles not of type I, D, or M. 
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and 



r]u 



K / 1 dr] vpscrr 2uar var^ 



(5.9) h = Q r^''^ 1^ 1 T^? - ^ + n \72 ^ .2, - ^ ) 

Jo r]\s\\l-ridx Tu [1 - ri){u^ + Q) u J 

/33 = Ci / -rr^^' C = id- 

Jo V\s\u 

Observe that a is periodic in e and that the /3j are ah real, with /3i and (3^ positive while the sign 
of (32 depends on the profile w. 

Remark 5.2. By Assumption 11.61 the stability function for the von Neumann step-shock, L,,7v(C)) 
is nonvanishing in 3ftC > 0. Stability functions for step-shocks were defined and analyzed in | Er4j . 
Remarkably, it turns out that 

(5.10) L,^(C) = Li(C), 

and the analysis of |Er4] shows that since Li is nonvanishing in > 0, in fact 

(5.11) Li(C) > for Ci > co??'/'(0+) 

and increases monotonically with Ci- The function in (jS.lUp is a nonvanishing multiple of the 
stability determinant for step-shocks defined in [M] , 



Theorem 5.3. Let t = C^ + t^, where |z^| < R. Under Assumptions \l.^ \1.3l and \1.4\ we have: 

a) For C £ {3^C ^ 0} \ Class III there exists a positive constant e{C,R) such that V{T,e) ^ for 
e > e{C,R)- This result holds even for profiles that are not of type I, D, or M . 

b) For a profile of type D and C G Class III+ \ {ico'q^^'^ (0+) , ic^o'q^^'^ (oo)} , there exists a positive 
constant e(C, R) such that V{t, e) / for e > e(C, ^). 

c) For a profile of type I, suppose that for iCi G Class ///+ \ {icQrf-^'^ {{)+), icQr]^/'^{oo)} one has 
the inequality 

(5.12) \L2{iCi)\^MP2{Ci)] < Li{iCi). 

Then there exists a positive constant e{Ci, R) such that V{iCie + i^,e) ^ for e > e{Ci, R)- 

d) For a profile of type I, suppose that for some id E Class ///+ \ {zco?7"^/^(0+), ico?/^''^(oo)} one 

has 

(5.13) \L2{iCi)\eMP2{Ci)\> Li{id). 

Let V satisfying {vl < R be an arbitrary point on the vertical line defined by 

(5.14) Li = |L2|exp[/32-K(z.)/33], 

and suppose now that R > for Ri^ as in ()5.14p . For any 5 > there exists a positive constant 
€{Ci, R, S) such that for periodically distributed e > e{Ci, R, S) satisfying 

e = [(2n-l)7r-9(i/)/33]//3i, L2 > 

(5.15) I 

e = [(2n + -)7r-9(i/)/33]//3i, L2 < 0, 
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V{iCi€ + I'je) has a zero within a distance 6 of v. 

(e) For a profile of type M, suppose thatiQ G Class ///+\{ico7?^''^(0+), zco??^''^(xAf ), «cor?^/^(oo)}. 

If there is a single turning point and d = -^3^(2;*) < 0, the conclusion of part (b) holds for 
C = iCi- 

Suppose there is a single turning point and d > 0. If inequality (|5.12p holds, then the 
conclusion of part (c) holds, while if inequality ()5.13p holds, then the conclusion of part (d) holds. 

Suppose there are two turning points X2, < xu. If inequality (15.120 holds, then the conclusion 
of part (c) holds, while if inequality ()5.13p holds, then the conclusion of part (d) holds. Recall that 
when there are two turning points, the contour in the integral defining a encloses only X2* ■ 

The constants e{C,R,5) (resp. e{Ci,R,5)) above can be chosen uniformly for (resp. Q) in 
compact subsets of the sets of values described in the respective cases (or subcases) listed above. 

Proof. 1. First observe that as a consequence of (15.40 . the estimates (15. 5p . (15. 6p hold with La 
replaced by L. Parts (a) and (b) then follow immediately from the nonvanishing of Li{C,) and parts 
(1) and (2) of Proposition ISTTl 

2. From (j5.8p we have 

(5.16) \a\ = exp(/32 - ^{v)p:i), /^s > 0. 

The inequality ()5.12p thus implies that Li^iQ) + a(e, iQ, v)L2{iC,i) does not vanish for "Slv > 0. Part 
(c) now follows from (15. 4p and part (3) of Proposition 15.11 

3. When the inequality (j5.13p holds, (|5.8p implies that the analytic function of v given by 

(5.17) V;(z.;0,e) := Li(iCi) + a(e, ^0, z^)^2«i) 

has zeros on the vertical line (|5.14p at the e values given by (|5.15p . The nonvanishing of (3^ and 
-^2(^0) implies these zeros are simple. Observe also that 

(5.18) |a^K(z^;Ci,e)l = e^2-»Hft^3|^2| is independent of e. 

Fix an arbitrarily small 5 > and let Cs be a circular contour centered at v of radius 5. For e large 
enough satisfying ()5.15p . part (3) of Proposition [5TT] and (|5.18p imply 

(5.19) \V{iCie + u,e)-Va{v;C^,e)\ < ^ < \Va{yXue)\ on Cs. 

Part (d) now follows from Rouche's Theorem. The subcases of part (e) follow in the same way from 
part (4) of Proposition 15.11 

4. The statements about uniform choices of e are direct consequences of the corresponding 
statements in Proposition 15.11 

□ 

Applications. 

1) We refer to |Er3] for a detailed discussion of the consequences of Theorem 15.31 For one 
reaction A ^ B detonations with the Arrhenius rate law (jl.30p . Erpenbeck first identifies ranges of 
various physical parameters for which profiles of type I, D, or M do occur. Then he uses Theorem 
15.31 to identify ranges of the parameters for which unstable zeros of V are actually present in cases 
I and D. Many of these results are also reported Chapter 6 of the book |FD] . 

2) Starting with an open interval of Q values satisfying (I5.13p . one obtains from (I5.15p . say 
when 9(1^) = 0, corresponding e intervals of unstable wavenumbers. For n sufficiently large these 
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intervals overlap and one finds that all wavenumbers above a certain cutoff (depending on the 
interval of C values) are unstable. 

3) Most of the conclusions of [Er3j for type M profiles do not rely on his analysis of the turning 
point problem for Cj,M := cor]^^'^{xM), where xm is the location of the maximum. In particular, by 
applying part (e) of Theorem 15.31 to values of Q near Ci,M, but not equal to it, Erpenbeck derives 
the algebraic condition 

T^/^ N I I drj vpsar 2uar var^] , , 

'-{—.l- — ^ (i-.)K + C|m) - — j " " 

as a sufficient condition for instability for profiles of type M. The inequality (j5.20p implies the 
existence of unstable zeros of V(iC^ie + z/, e) for C^i near Ci,M- A separate criterion is derived in |Er3j . 
p. 96, for V to have unstable zeros exactly at the special value Ci,M@rhis condition cannot be 
derived from Theorem! 



A Appendix 

A.l A variable-coefficient gap lemma 

In this section we state and prove Lemma lA.H which is our main tool for determining the 
asymptotic behavior as e = ^ — )■ oo of the decaying solution 9{x,T,e) of ()1.7p on [a, oo). This 
lemma was introduced in [Z2] and used there for the study of stability of ZND detonations with 
respect to high frequency one-dimensional perturbations. 

Consider a first-order system 

(A.l) V' = AP{x; h)V := MP{x; h)V + 6^(2;; h)V on x > 0, 

with V{x) G C^, depending on a parameter p ^ P C and a distinguished small parameter 
h>0. 

Lemma A.l ( |Z2] ). Suppose there exist positive constants /3, C such that for all x > and p G P: 
(A.2) \eP{x,h)\<Ch^e-l^^'' and 

(A.3) mPix, h) > -{h5P{h) + C/ie"^^^), where < 5P{h) < (5, < /3 

for all /il^ Assume further that there exists a nonzero constant vector such that 
(A.4) MP{x, h)VP = for all x>0. 

Then there exists an ho independent of p € P such that for < h < ho, there exists a solution 
VP{x, h) of (lA.ip defined on x > satisfying 

(A.5) \{VP{x) - VP)\ < Cihe~^*''''\VP\ for x>0. 



For ^i,M a formula like (|4.33|) is derived in [Er3| . but with a completely different a. 
Here KM := i(M + M*). 
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Proof. 1. We seek, equivalently, a solution of the integral fixed-point equation 

(A.6) TVix) = Vr + r ^y^^QPiy, h)V{y)dy, 

J +00 

where J^y^^ is the solution operator of V = M^V from y to x; that is, W{x) := F^~^^w{y) satisfies 
(A.7) W^ = MPW, WU=y = w{y). 

2. We claim 

(A.8) 11-^^^"^ II < Ce^^^'-^'^^y-'''^ < Ce''^*^y-^^ for y > x. 

Indeed, for W{x) as in (\A.7\i we have using ()A.3P 

d^(|l^|2) = 2^{W,W^) = 2{W,?R:MPW) > -2(MP(/i) + C/ie-^^^)|T^|2 ^ 
^"^■^^ |WU/|Ty| > -(MP (/i) + C/ie"'^'^^), 

and integrating this inequality from x to y gives 

POO 

(A.IO) |VF(x)|/|u;(y)| < C7ie'*'''('^)(^-^), where Ci = / Che-^^^'dx. 

Jo 

3. For /i > sufficiently small, this implies that 7" is a contraction on L°°[0, 00). For, applying 
()A.2p and (jA.SP we have (with |y|oo = |^|l°°[o,oo)) 



;>oo 

(A.ll) \rVi - TV2\ (x) < Ch^\Vi - V2\oo / e'^*'^(^-")e-^'^^dy < C"/i|yi - ^2! 



— (5,/ix 



which for h sufficiently small is less than ^|Vi — V2|oo- By iteration, we thus obtain a solution 
VP G L°°[0, 00) of y = TV. Further, taking Vi = VP, ^2 = in (|XTT]l we obtain both 

\VP-VP\{x)<C'h\VP\^e-'''''', 

which together imply (|A.5p . □ 

Remark A. 2. Strictly speaking. Lemma lA.ll is a "weak" gap lemma realizing only the fastest 
decaying mode of ()A.ip . whereas the full gap lemma construction of |GZ| (see, for example. Cor. 
2.4, [GZ], or Prop. 3.1 of |ZHj ) yields intermediate modes as well. However, a straightforward 
modification of (]A.6P following that of Prop. 3.1, [ZH] yields intermediate modes as well, provided 
that MP splits into diagonal blocks such that 3f?M^(x,/i) ^ -{h5P{h) + Che'l^^''). 

Remark A. 3. On an interval [a/h,b/h], corresponding in coordinates x = hx to a finite interval 
[a, 6], essentially the same argument as in the proof of Lemma I A . 1 1 applies with modified mapping 

(A.13) TV{x) = VP + r J^y^^QPiy, h)V{y)dy 

Jb/h 

to yield an analogous result 1(^^(2;) — Kf)| < Ci^|Kif| for x G [a/h,b/h], assuming only uniform 
boundedness of the coefficient matrix AP, \@P[x,h)\ < Ch'^, and ^MP[x,h) > —Ch, or, for inter- 
mediate modes, ^M^{x,h) ^ —Ch, similarly as discussed in Remark |A.2[ Rescaling x — )• hx, and 
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comparing with Theorem 13. ll fb). we see that, restricted to a finite closed real interval, the Method 
of Parameters and the variable-coefficient gap lemma constructions thus give essentially the same 
results (here restricted to level m = 0). On an infinite interval, on the other hand, it is readily 
checked that the truncation error |$m — ^\ resulting from the Method of Parameters construction 
described in (jl.lSp generically does not decay in x, even for e<I>(e, z) = e^Q + <I>i with <I>j constant 
(rather than just asymptotically constant), and so the method of parameters in general fails 1^ 

A. 2 ZND equations 

The ZND equations for the unknowns {v, u, S, A) (specific volume, particle velocity u = {ux,Uy, Uz), 
entropy, and mass fraction of reactant) are 

dtv + u • V?; — vV • u = 
dfU + u • Vu + vVp = 
^ ■ ' dtS + u-VS = -rAF/T:=<i> 

dtX + u-VX = r, 

where p = p{v, S, A) is pressure, T is temperature, AF is the free energy increment, and r{v, S, A) 
is the reaction rate function. 

A. 3 The stability function 1/(t, e). 

The stability function defined in |Erlj and used in [Er3] is given by 

(A.15) V{t, e) = T6i(r, e) + ieb2{T, e) - 6(0, r, e) • (rht + iehy), 

where (with v' = ^ and m = u/v, a constant independent of x) 

(A.16) 

/•oo 

gt = -{v',u',0,S',X'Y, gy = (0,0, -vp',0,OY, b, = - eix,T,e) ■ A-\x)gjix)dx for j = 1,2, 

Jo 



(A.17) 



( 2{l -'n+)g+/m \ 
T+7?+ + 2(1 - 


—m{v^ — 




and hy the single has nonzero component {hy)^ = m{v-—v^). Here v±, for example, are components 
of the states 'w± := 'w{0±) just to the right and left of the von Neumann shock, and 



(A.18) 



9+ 



In section [5^ we work with the simpler form of V (j5.3p in which the integals bj do not appear. 

The functions -Z^i(C) and L2{C), appearing in Proposition 15.11 and used to compute zeros of V, 
are given explicitly in [Er3] as: 



Specifically, the righthand side of p.l2p blows up as z — s- oo. 
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(A.19) 



^i(C) 



^-(1 - Xv) 



^+C(C + ^+g+) 



^ = 2 - (1 - r]){l - Xv)v-Ps/T, Xv = v+/v^, 
where ± denotes evaluation at 0±. 



(A.20) 



^2(0 



UiU- 



U+U- 



A.4 CoefRcients appearing in the linearized systems 

The matrix coefficients appearing in the reduced system (jl.Sp are 



A, 



(A.21) 



/ u 


-V 





\ 




1 





— V 
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VPv 
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vps 


vpx 






















u 








, Ay — 


VPv 








VPS 


vp\ 








u 






















V 








u J 




V 











/ 



B 



V 



v{ci/v^y 





V 



V 

u' 

S' 
X' 













vp's 


-rs 



\ 

vp'^ 





where (') denotes differentiation with respect to x and Cq = —v'^pv{v,S,X). These matrices are 
obtained from the corresponding matrices in the unreduced 6x6 system by deleting the fourth row 
and fourth column of each of the latter matrices. In each of the unreduced matrices Ay (coefficient 
of a) and B, the fourth row and fourth column consist only of zeros, while the fourth row and fourth 
column of Ax have the fourth component u and all other components 0. The matrix coefficient of 
/? in the unreduced system is 



(A.22) 



A.. 



1 ° 








— V 





\ 






































VPv 











vps 


VPv 




















V 














) 



The matrix ^q{x,C,) in the transposed system (|1.9p is computed in |Er3] . p. 112 to be 



(A.23) 



^o(x,C) 



/ (i-'?)C 
™ 

(l-r?)C 
rjmu 

ii-v)ps( 



(1-'?)PaC 



m( 

rm 

(i-^)C 



V 

{l-ri)psC 

rjmu 
(1-'?)PaC 

rjmu 



im 
l-r) 





o\ 











c 

u 








¥s 
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u 


m 
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The matrix E = T ^<l>iT — T (|1.14p has En component given by 

n / A N — 'n)vps 1— "^/^ N 2 — 71 dri 

kC f 1 dr] VPS I/O , ^^ , ^^C + s C Idlns 
err [2i' + var^) 1 H err 



2rys \1 — 7] dx Tu u J ( + ks rjus 2 dx 

where a is defined ( |FD| .p.95) by 

(A. 25) a = v{dp/dX)e^^/cQ (e is specific internal energy). 

The formula for E22 is the same as (IA.24p . except that s is replaced by —s. 
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